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ABSTRACT 


DEANE  is  a  special-purpose  computer  language  designed  for  use 
in  a  time-shared  environment  by  a  ballistic  missile  defense  sys¬ 
tems  analyst.  In  essence,  it  is  a  sophisticated  calculator  whose 
modular  design  allows  the  user  to  request  basic  computations  in 
a  convenient  fashion.  The  interpretation  of  the  computational  re¬ 
sults  is  not  made  by  DEANE  under  the  pretense  of  being  a  sim¬ 
ulator,  but  is  left  to  the  user. 

This  report  is  a  user’s  manual  describing  the  computations  avail¬ 
able  and  giving  examples  of  how  they  may  be  ordered  to  solve 
typical  problems.  Also  presented  is  a  description  of  the  logical 
and/or  mathematical  foundations  of  the  computations. 


Accepted  for  the  Air  Force 

Joseph  R.  Waterman,  Lt.Col.,  USAF 

Chief,  Lincoln  Laboratory  Project  Office 


CONTENTS 


Abstract  iii 

I.  Introduction  1 

II.  System  Overview  2 

III.  Trajectory  Specification  Commands  4 

IV.  Trajectory  Observation  Commands  11 

V.  Trajectory  Interception  Commands  19 

VI.  "Cookie- Cutting”  Commands  29 

VII.  Convenience  Commands  3  3 

VIII.  Command  Summary  36 

IX.  Mathematics  and  Logic  38 

A.  Basic  Coordinate  System  38 

B.  KEPLER  38 

C.  EJECT  50 

D.  R_I_AXIS  52 

E.  INTEGRATE  55 

F.  PRINT  56 

G.  VISIBLE  61 

H.  INTERCEPT  64 

I.  CONTOUR  64 


DEANE:  A  COMPUTER  AID 
FOR  BALLISTIC  MISSILE  DEFENSE  ANALYSIS 


I.  INTRODUCTION 

DEANE  is  a  special-purpose  computer  language  designed  to  aid  in  certain  aspects  of  ballistic 
missile  defense  systems  analysis.  It  is  designed  to  help  solve  problems  which  involve  under¬ 
standing  primarily  geometric  and  kinematic  relationships.  DEANE  is  intended  to  operate  in  a 
time-shared  environment,  but  it  can  also  be  executed  under  a  batch  monitor. 

The  earth  is  modeled  as  a  sphere  which  may  be  either  stationary  or  rotating.  Trajectories 
are  initially  calculated  to  be  Keplerian;  however,  there  is  an  option  for  replacing  the  re-entry 
part  of  the  trajectory  with  values  attained  by  integration  of  the  equations  of  motion  using  a  real 
atmosphere  and  a  finite  ballistic  coefficient  for  the  re-entry  object.  These  Keplerian  trajecto¬ 
ries  may  be  determined  from  either  launch  site  and  impact  site  specification  or  by  specification 
of  ejection  velocity  components  from  a  base  trajectory.  Axisymmetric  objects  may  be  given  a 
cross  section  as  a  function  of  aspect  angle  by  specification  of  a  table,  and  their  orientation  may 
be  modeled  to  correspond  to  either  constant  rate  inplane  tumbling  or  a  spin-stabilized  situation 
with  coning. 

With  the  exception  of  the  parameters  necessary  to  calculate  the  signal -to-noise  ratio,  each 
sensor  is  modeled  only  for  geometric  purposes  by  its  position,  azimuth  boresight,  and  a  limit¬ 
ing  range,  elevation,  and  azimuth.  Each  interceptor  farm  is  modeled  by  its  position  and  what 
type  of  interceptor  is  stationed  there,  the  interceptor  type  referring  to  a  particular  set  of  fly¬ 
out  characteristics  as  specified  by  a  table  of  flyout  times  vs  tangent  range  and  altitude. 

Within  the  model  defined  above,  DEANE  can  perform  various  computations.  These  compu¬ 
tations  are  just  that,  leaving  any  interpretation  of  the  results  to  be  made  by  the  user,  not  by 
DEANE.  Typical  examples  would  be.  (a)  calculating  ideal  observables  on  a  trajectory  over 
some  time  interval  from  specified  sensors,  (b)  calculating  when  an  object  is  first  visible  to  a 
specific  sensor  as  determined  by  the  trajectory  and  a  limiting  range,  azimuth,  and  elevation 
associated  with  that  sensor,  (c)  calculating  the  earliest  possible  intercept  and  the  intercept  with 
the  latest  launch  from  some  set  of  farms,  subject  to  various  restrictions  such  as  a  minimum  in¬ 
tercept  altitude,  a  maximum  commit  altitude,  and  the  requirement  that  the  launch  must  occur 
at  least  so  many  seconds  after  any  one  of  a  set  of  specified  radars  has  been  able  to  detect  the 
object,  and  (d)  determining  a  ground  footprint  of  interceptor  coverage  subject  to  restrictions 
similar  to  those  above. 

Also  available  are  "cookie-cutting"  functions  for  investigating  the  deployment  of  sensors, 
farms,  and  Standard  Metropolitan  Statistical  Areas  (SMSA’s).  Typical  calculations  would  be  to 
determine  which  SMSA's  are  within  a  certain  distance  d^  of  any  farm  and  within  a  certain  dis¬ 
tance  d  of  any  sensor,  or  to  determine  which  farms  or  sensors  are  within  a  certain  distance 
of  specified  SMSA's. 

As  DEANE  is  intended  for  use  in  a  time-shared  environment,  significant  consideration  was 
given  to  making  the  specification  of  program  options  as  convenient  for  the  user  as  possible.  The 
resulting  solution  was  that  all  the  user  input  to  DEANE  is  in  the  form  of  a  command,  or  key¬ 
word,  followed  by  the  necessary  input  data  needed  by  that  particular  command.  Accordingly, 
DEANE  can  be  considered  as  either  a  program  or  a  special-purpose  language,  the  latter  being 
preferred  for  purposes  of  explanation,  understanding,  and  usage. 
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II.  SYSTEM  Q\  ERVTEW 


An  oversimplified  view  of  a  computer  will  help  understand  how  DEANE  works.  We  may 
consider  the  computer  to  be  some  organized  collection  of  memory  cells  and  a  set  of  hardware 
that  can  perform  various  arithmetic  and  logical  operations  on  these  cells.  Each  cell  is  made 
accessible  to  the  arithmetic/logical  hardware  by  specifying  a  bit  pattern  address  which  at  the 
user  level  is  equivalent  to  specifying  a  mnemonic.  Whatever  is  in  one  of  these  cells  is  unchanged 
each  time  the  contents  are  read  for  computational  purposes,  and  the  contents  are  changed  only 
when  something  else  is  written  into  the  cell.  Accordingly,  there  are  two  kinds  of  commands  in 
DEANE:  (a)  data  defining  commands  which  cause  a  sequence  of  numbers  to  be  read  and  stored 
in  appropriate  cells,  overwriting  the  previous  contents;  and  (b)  action  commands  which  initiate 
the  desired  calculations  using  the  present  values  resident  in  the  appropriate  cells.  Data  defin¬ 
ing  commands  are  the  only  means  by  which  certain  values  (e.g.,  trajectory  launch  site)  can  be 
changed.  Action  commands  either  produce  output  or  determine  values  for  variables  which  need 
be  known  only  to  the  computer  (e.g.,  the  argument  of  the  trajectory  apogee). 

A  statement  in  DEANE  consists  of  a  command  followed  by  a  sequence  of  numbers  and/or 
mnemonics.  This  data  sequence  varies  in  length  and  meaning  for  the  different  commands.  The 
leading  keyword  of  each  statement  must  not  be  in  quotes  and,  except  where  noted,  the  entire 
statement  must  be  entered  on  one  7  2-space  line.  Either  one  comma  and/or  any  number  of  blanks 
may  be  used  to  separate  the  entries  within  a  statement.  The  units  for  the  input  for  all  the  com¬ 
mands  are  identical;  time  is  in  seconds,  distance  is  in  kilometers,  speed  is  in  kilometers/sec¬ 
ond,  angles  are  in  degrees,  angular  rates  are  in  degrees/second,  latitude  is  in  degrees  north, 
and  longitude  is  in  degrees  east. 

To  help  clarify  the  above  comments  and  to  give  insight  into  how  computations  can  be  achieved 
using  the  entire  instruction  set,  Fig.  1M  is  an  example  of  some  of  the  commands  and  how  they 
might  be  used. 

1  -3-13045  I 

launch  38.0  110.0  /  ,  foe 

Impact  40.8,-75.0  5.  friend 
parameter  0  1 

kep 1 e  r 

TOTAL  FLIGHT  TIME  IS  2070.48  SECONDS, 
a  1 1 1 tude  70.0 

THE  OBJECT  IS  AT  ALTITUDE  70.00  AT  TIME  -26.43. 

parameter  1,2,10 
kepi er 

TOTAL  FLIGHT  TIME  IS  1710.14  SECONDS, 
a  1 t i tude  70 

THE  OBJECT  IS  AT  ALTITUDE  70.00  AT  TIME  -51.19. 

stop 

T=3.01/6.96  13:52:58 

Fig.  11-1.  Example  of  things  to  come. 

The  first  statement  fixes  the  site  named  'FOE',  located  at  38.0°  north  latitude,  110.0°  east 
longitude,  and  0  kilometers  height  as  the  launch  site.  The  second  statement  specifies  that  im¬ 
pact  is  to  be  at  5.0  kilometers  over  40.8°  north  latitude  and  —75.0°  east  longitude,  and  this  site 
is  given  the  name  'FRIEND'.  Statement  three  indicates  that  a  minimum  energy  trajectory  is 
desired  for  a  nonrotating  earth,  and  the  statement  'KEPLER'  then  causes  calculation  of  the  pa¬ 
rameters  specifying  that  trajectory  and  the  line  of  output  that  follows.  The  next  command  pro¬ 
duces  the  line  of  output  indicating  that  the  object  is  at  70.0  km  altitude  26.43  seconds  before  im¬ 
pact.  'PARAMETER  1,  2,  10'  sets  the  parameters  for  future  trajectory  definitions  to  be  a 
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rotating  earth  with  a  specified  impact  velocity  angle  of  10°.  ’KEPLER’  then  determines  the 
trajectory  using  the  launch  site  ’FOE’,  the  impact  site  ’FRIEND’,  a  rotating  earth,  and  an  im¬ 
pact  velocity  angle  of  10°.  ’ALTITUDE  70.0'  then  produces  the  following  line  of  output,  and 
’STOP'  terminates  the  execution. 

In  this  example  the  commands  ’LAUNCH',  'IMPACT',  and  'PARAMETER'  are  data-defining 
commands.  ’KEPLER’  is  an  action  command  which  determines  a  set  of  values  which  need  be 
known  only  to  the  computer  and  produces  one  line  of  output.  ’ALTITUDE’  is  an  action  command 
which  serves  only  to  produce  output. 

The  complete  list  of  instructions  is  presented  in  Secs. Ill  through  VII.  The  description  of 
each  command  is  separated  into  three  sections:  USE,  which  describes  why  the  command  is  in¬ 
cluded  in  DEANE  and  how  it  might  be  used,  RESPONSE,  which  describes  the  computer's  re¬ 
sponse  to  the  command,  and  ERROR,  which  describes  syntactic  and/or  semantic  limitations  on 
the  command. 

For  all  the  commands,  the  error  message  'ILLEGAL  INPUT  INVOLVING  THE  SYMBOL 
"x"  '  will  be  printed  when  DEANE’s  input  routine  recognizes  an  illegal  representation  for  a 
number  in  a  statement.  The  error  message  'SOME  INPUT  VALUE  IS  UNACCEPTABLE. 
PLEASE  RETYPE  THIS  LINE.'  will  be  printed  when  the  value  of  some  input  is  syntactically  un¬ 
acceptable.  Usually,  all  input  up  to  the  point  of  infraction  is  accepted,  and  all  input  including 
and  past  the  illegal  number  is  not  accepted.  The  user  response  should  be  to  retype  the  entire 
command  in  which  the  error  occurred. 
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III.  TRAJECTORY  SPECIFICATION  COMMANDS 


LAUNCH  lat  long  ht  mnemonic 

USE  -  This  specifies  the  launch  site  of  a  trajectory. 

RESPONSE  -  The  launch  site  is  set  to  flat’  degrees  north  latitude,  'long*  degrees  east 
longitude,  and  'ht*  kilometers  above  the  earth.  Any  'mnemonic'  entered  after  'ht'  is 
ignored  by  DEANE,  but  it  may  serve  the  user  for  purposes  of  identification.  If  'ht'  and 
'mnemonic1  are  omitted,  'ht'  is  set  to  zero. 

ERROR  -  None. 

IMPACT  lat  long  ht  mnemonic 

USE  -  This  specifies  the  impact  site  of  the  Keplerian  trajectory. 

RESPONSE  -  The  impact  site  is  set  to  'lat'  degrees  north  latitude,  'long'  degrees  east 
longitude,  and  'ht'  kilometers  above  the  earth.  Any  'mnemonic'  entered  after  'ht'  is 
ignored  by  DEANE,  but  it  may  serve  the  user  for  purposes  of  identification.  If  'ht'  and 
'mnemonic'  are  omitted,  'ht'  is  set  to  zero. 

ERROR  -  None. 

PARAMETER  rotate  type  valuel  value2 

USE  -  This  command  specifies  the  information  additional  to  the  launch  and  impact  sites 
necessary  to  allow  calculation  of  a  Keplerian  trajectory. 

RESPONSE  -  If  'rotate'  equals  one,  the  earth  is  modeled  as  rotating;  otherwise  it  is 
modeled  as  fixed.  There  are  four  different  types  of  parameters  which  are  indexed  by 
'type'  and  given  values  via  'valuel'.  The  association  is  given  in  Table  III  — 1. 


TABLE  III-l 

INTERPRETATION  OF  DATA  FOR  COMMAND  'PARAMETER' 

'type' 

Interpretation  of  'valuel' 

1 

None  —  minimum  energy  calculation. 

2 

Angle  between  local  horizontal  and  velocity 
vector  at  impact. 

3 

Velocity  magnitude  at  impact.  Since  this  leads 
to  two  possible  trajectories,  a  positive  value 
yields  the  lofted  trajectory  choice  and  a  negative 
value  yields  the  depressed  trajectory  choice. 

4 

Apogee  height  in  kilometers. 

If  the  trajectory  specified  is  an  orbit,  'value2'  represents  the  number  of  complete 
orbits  desired  before  satisfying  the  impact  conditions. 

ERROR  -  'type'  must  be  an  integer  between  1  and  4,  and  'rotate'  and  'value2'  must  be 
integers. 
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KEPLER  and  KEPLERP 


USE  —  At  any  time,  the  user  has  access  to  only  one  trajectory  for  purposes  of  per¬ 
forming  computations.  This  command  provides  one  means  of  attaining  such  a  working 
trajectory. 

RESPONSE  -  The  Keplerian  trajectory  as  specified  via  the  most  recent  values  from 
’LAUNCH',  ’IMPACT’,  and  'PARAMETER'  is  determined.  'KEPLERP'  will  cause 
printing  of  various  characteristics  about  the  trajectory  (longitude  of  the  ascending  node, 
inclination,  argument  of  apogee,  eccentricity,  semi-latus  rectum,  semi-major  axis, 
flight  ground  range,  total  time  of  flight,  and  the  velocity,  azimuth,  and  flight  path 
angle  at  launch  and  impact).  For  'KEPLER',  only  the  total  time  of  flight  is  printed. 

At  the  successful  completion  of  either  of  these  commands,  the  working  trajectory  is 
the  one  calculated,  overwriting  any  previous  working  trajectory. 

ERROR  -  There  are  several  situations  where  the  data  specified  via  'PARAMETER' 
are  not  compatible  with  the  launch  and  impact  sites.  If  the  velocity  impact  angle  is  too 
large,  a  message  will  be  printed  and  a  trajectory  with  an  apogee  of  25,000  km  is  sub¬ 
stituted.  If  the  impact  velocity  is  too  small,  a  trajectory  with  minimum  launch  velocity 
is  substituted  and  a  note  is  printed  to  that  effect.  Should  the  impact  velocity  be  too 
large,  a  message  is  printed  and  an  appropriate  substitution  is  made  for  the  impact 
velocity.  If  computational  problems  arise,  the  trajectory  calculation  is  abnormally 
terminated  and  a  message  to  that  effect  is  printed. 

Example  1 

Of  recurring  interest  in  many  analyses  is  the  relation  of  re-entry  angle  to  ground  range. 
Using  only  the  DEANE  commands  introduced  thus  far,  we  can  determine  the  re-entry  angle  as 
a  function  of  trajectory  ground  range  for  a  given  burnout  velocity.  A  sequence  of  commands  to 
do  this  might  be: 

impact  48.0  —80.0  0.0 
parameter  0  3  7.2 
launch  70.0  120.0  200.0 
keplerp 

launch  67.0  120.0  200.0 
keplerp 

launch  64.0  120.0  200.0 
keplerp 


launch  19.0  120.0  200.0 
keplerp 

parameter  0  3  —7.2 
launch  70.0  120.0  200.0 
keplerp 

launch  67.0  120.0  200.0 
keplerp 


The  launch  velocity,  ground  range,  and  impact  angle  printed  as  a  result  of  'KEPLERP'  then 
specify  the  relationship  of  interest.  (By  conservation  of  energy,  constant  impact  altitudes  and 
velocities  with  constant  launch  altitudes  give  constant  launch  velocities.)  Figure  III- 1  shows 
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Fig.  Ill  - 1 .  Re-entry  angle  vs  ground  range  as  calculated  by  DEANE. 

such  curves  for  several  burnout  velocities,  from  which  the  analyst  could  make  observations 
about  possible  re-entry  angles  for  a  particular  booster  and  ground  range  or  about  the  effects  of 
changing  the  payload  weight  and  burnout  velocity  and  hence  range  of  re-entry  angles. 


EJECT  time,  delvx,  delvy,  delvz,  alt 

USE  -  This  command  provides  another  means  of  obtaining  a  working  trajectory  to 
model  decoy  ejection,  tank  displacement,  or  deorbiting. 

RESPONSE  -  At  the  time  'time1  before  impact,  the  present  working  trajectory's  veloc¬ 
ity  vectors  are  incremented  by  'delvx',  'delvy',  and  'delvz'  km/sec  and  calculations 
then  determine  a  new  Keplerian  trajectory.  The  impact  point,  velocity,  angle,  and 
time  to  impact  of  the  resulting  trajectory  are  printed,  where  impact  is  taken  to  be  at 
the  altitude  'alt'.  The  coordinate  system  used  for  specifying  the  velocity  increments 
is  seen  in  Fig.  Ill- 2;  x  is  the  inplane  tangent  to  the  old  trajectory  at  the  time  'time', 
y  is  the  inplane  normal  to  x,  and  z  is  the  right-hand  normal  to  x  and  y. 


Fig.  Ill  -  2.  Coordinate  system 
for  input  to  'EJECT'. 


ERROR  -  If  the  new  trajectory  never  reaches  the  altitude  'alt',  an  arbitrary  acceptable 
altitude  for  impact  is  chosen,  and  a  message  is  printed  to  this  effect. 
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R_I_AXIS  time,  lat,  long,  alt 

USE  -  This  command  aids  the  user  in  selecting  input  values  for  the  command  'EJECT*. 

RESPONSE  —  The  vector  (in  the  coordinate  system  sketched  in  Fig.  III-2)  along  which 
small  velocity  increments  at  time  'time*  will  produce  impact  at  'lat*  latitude,  'long* 
longitude,  at  a  height  'alt'  over  the  earth  is  determined  and  printed.  Also  printed  are 
the  magnitude  and  components  of  the  minimum  magnitude  velocity  increment  necessary 
to  achieve  this  displacement.  The  calculations  are  based  on  linear  approximations  and 
a  nonrotating  earth.  They  are  thus  valid  only  for  desired  impact  points  near  that  of 
the  present  trajectory. 

ERROR  -  Because  of  the  nature  of  the  computation,  if  the  working  trajectory  never 
reaches  the  altitude  'alt',  the  command  is  terminated. 


Example  2 


Suppose  we  desire  to  model  a  decoy  whose  impact  is  0.5°  (about  55  km)  uprange  from  the 
impact  point  for  a  base  trajectory.  The  sequence  of  DEANE  commands  with  their  output  (seen 
in  Fig.  Ill  -  3 )  would  achieve  this.  The  first  four  commands  specify  and  cause  calculation  of  a 
base  trajectory.  The  *R_I_AXIS'  command  indicates  the  displacement  can  be  achieved  at  time 
—  1 840  with  a  velocity  increment  along  the  vector  (—.01 0,  .  022,  0)+c(— .  207,  2.  86,  0)  where  c  can 
be  arbitrarily  chosen.  Choosing  the  minimum  velocity  magnitude  from  *R  I_ AXIS'  as  input  for 
'EJECT*  causes  calculation  of  this  trajectory  which  does  indeed  impact  0.5°  uprange  of  the  base 
trajectory. 

|  -3  -  1304  8  | 

launch  42.1  83.0  200.0 
I mpac t  48.0  -97.0  0.0 
parameter  0  2  20.0 
kep 1 e  r 

TOTAL  FLIGHT  TIME  IS  1859.50  SECONDS. 
r_i_axls  -1840  48.5  -97.0  0.0 

DELVX  =  -1.010E-02  +  -2.070E-01  TIMES  C. 

DELVY  =  2.215E-02  +  2.865E  00  TIMES  C. 

DELVZ  =  6.481E-06  +  2.821E-07  TIMES  C. 

FOR  MINIMUM  VELOCITY  MAGNITUDE  OF  8.474E-03 

DELVX  =  -8.452E-03,  DELVY  =  -6.107E-04,  DELVZ  *  6.478E-06 

eject  -1840  -8.452e-3  -6.107e-4  6.478e-6 

IMPACT  AT  -0.00  KM  ABOVE  48.501  DEG  LAT  AND  -97.000  DEG  LONG 
WITH  A  SPEED  OF  7.2771  KM/SEC  AND  AT  AN  ANGLE  OF  20.01  DEGREES. 

FLIGHT  TIME  FROM  VELOCITY  INCREMENT  IS  1829.37. 


Fig.  Ill- 3.  DEANE  commands  for  Example  2. 

To  extend  this  example,  the  minimum  velocity  magnitude  and  associated  change  in  flight 
time  can  be  calculated  as  a  function  of  the  inplane  displacement  by  repeating  the  'R_I_AXIS'  and 
'EJECT'  sequence  for  several  displacements.  The  results  of  doing  this  for  a  10,000-km  ground 
range  trajectory  are  seen  in  Fig.  M-4  where,  for  all  three  re-entry  cases,  ejection  is  about 
20  sec  after  200-km  altitude. 
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INPLANE  DISPLACEMENT  (km) 


Fig.  Ill -4.  Minimum  velocity  increment  and  flight  time  increase  vs  inplane  displacement 
as  calculated  by  repetition  of  ,R_I_AXIS’  and  'EJECT*. 


BETA_TABLE  / 3 v  hy  py  h,,  py  hy  P4 

USE  -  This  command  defines  a  ballistic  coefficient  (p)  history  as  a  function  of  altitude. 
This  is  necessary  input  for  achieving  trajectories  whose  re-entry  portion  is  the  result 
of  integration. 

RESPONSE  -  The  approximation  to  the  p  history  is  defined  by  the  piecewise  step  function 
shown  in  Fig.  III-5.  The  /T  units  are  taken  as  pounds  per  square  foot.  If  any  *h.f 
pairs  are  omitted  in  the  data  sequence,  'In'  is  assumed  to  be  zero. 

ERROR  -  fh.  ,  , '  must  be  less  than  'h.'. 

1 4 1  1 

“  18-3-13050 


^ ^ 

Fig.  Ill  -  5 .  Interpretation  of  input 
data  for  'BETA  TABLE'. 


— - INCREASING  'h  '  'h  '  'h  '  0 

1  2  3 

ALTITUDE  (km) 


'H 


■a; 


£(p*f) 


INTEGRATE  alt,,  alt, 

1  2 

USE  -  This  is  the  third  means  by  which  working  trajectories  are  attained,  and  it  is 
intended  to  provide  better  modeling  of  the  re-entry  portion  of  a  trajectory. 

RESPONSE  -  The  re-entry  part  of  the  present  working  trajectory  is  replaced  by  values 
attained  from  integrating  the  equations  of  motion  using  a  real  (tabular)  atmosphere  and 
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the  ballistic  coefficient  history  specified  via  'BETA_TABLE'.  The  integration  starts 
at  'alt/  km  altitude  and  terminates  at  'alt/  km  altitude  with  the  initial  state  vector 
determined  from  the  Keplerian  specification  of  the  working  trajectory.  As  drag  can 
significantly  alter  the  impact  point  of  such  a  trajectory  from  the  vacuum  impact,  the 
new  impact  point,  speed  at  impact,  and  impact  velocity  angle  with  the  local  horizontal 
are  printed.  Also  printed  is  the  time  increase  between  ' alt^ '  and  1  alt ^ 1  km  altitude  due 
to  the  integrated  trajectory. 

ERROR  -  If  the  integrated  trajectory's  time  to  ' alt ^ T  km  from  ' alt ^ 1  km  altitude  is 
greater  than  ZOO  sec,  the  integration  is  terminated  before  completion  and  a  message  is 
printed  to  that  effect. 

Example  3 

The  slowdown  time  is  often  of  interest  in  analyses  and  can  be  calculated  using  the  two  com¬ 
mands  presented  above.  Figure  Ill  -  6  shows  a  set  of  DEANE  commands  and  their  output  which 
indicate  that  for  a  10,000-km  ground  range  trajectory  and  a  constant  /3  of  600  there  are  9.6  sec 
of  slowdown  between  100-  and  3-km  altitude  for  a  30°  re-entry  trajectory.  For  a  10°  re-entry 
trajectory,  there  are  43.7  sec  of  slowdown.  By  executing  'INTEGRATE'  with  several  different 
constant  /3's  and  three  different  re-entry  angles,  the  relationships  shown  in  Fig.  Ill -7  may  be 
determined. 
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IV.  TRAJECTORY  OBSERVATION  COMMANDS 


PRINT  <p±  <p2  <?3..  .  <p  END 

USE  —  This  command  allows  the  user  to  specify  which  sensor  observables  are  to  be 
printed. 

RESPONSE -The  mnemonics  '<Py  <p^9  (py  ...  tp  '  specify  what  radar  observables  are 
to  be  printed  when  the  request  for  these  values  is  made.  This  command  may  extend 
over  more  than  one  line  and  is  terminated  by  the  mnemonic  'END'.  The  list  of  mne¬ 
monics  and  their  interpretations  is  given  in  Table  IV- 1,  and  the  default  set  for  this 
command  is  the  top  11  in  this  table. 


TABLE  IV-1 

INTERPRETATION  OF  INPUT  MNEMONICS  FOR  'PRINT' 

Mnemonic 

Interpretation 

TIME 

time  in  seconds 

RANGE 

sensor  range  in  kilometers 

RDOT 

sensor  range  rate  in  km/sec 

RDDOT 

sensor  range  acceleration  in  km/sec^ 

ELE 

elevation  in  degrees 

ELEDOT 

elevation  rate  in  deg/sec 

AZIMUTH 

azimuth  in  degrees 

AZIDOT 

azimuth  rate  in  deg/sec 

ALTITUDE 

altitude 

GRANGE 

ground  range  from  sensor  in  kilometers 

V A ANGLE 

velocity  aspect  angle  in  degrees 

BAANGLE 

body  aspect  angle  in  degrees 

SIGMA 

body  cross  section  in  square  meters 

S/N 

signal/noise  ratio  in  decibels 

LAT 

latitude  of  object  in  degrees  north 

LONG 

longitude  of  object  in  degrees  east 

ERROR  —  If  any  mnemonic  entered  is  not  in  Table  IV- 1,  a  message  to  that  effect  will 
be  printed  and  that  mnemonic  will  be  ignored.  If  too  many  mnemonics  are  listed  (as 
limited  by  132  characters  for  the  output)  a  message  will  be  printed  and  the  command 
will  be  terminated. 

POINTING  lat,  long,  alt,  rate,  angle 

USE  —  This  command  sets  up  the  parameter  values  for  modeling  spin  stabilized  orien¬ 
tation  with  coning  to  allow  computation  of  the  body  aspect  angle. 

RESPONSE  —  A  vector  in  inertial  space  pointing  toward  Tat’  degrees  north  latitude  and 
'long'  degrees  east  longitude  on  the  earth's  surface  from  the  trajectory  point  specified 
by  'alt'  is  established  as  the  axis  of  orientation.  If  'altr  is  positive,  the  trajectory  point 
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is  'alt'  km  altitude  on  the  impact  side  of  apogee;  if  it  is  negative,  the  point  is  J 1  alt 1  | 
km  altitude  on  the  launch  side  of  apogee.  The  object  is  then  given  a  coning  rate  of 
'rate'  dcg/sec  at  an  angle  of  'angle'  degrees  with  the  axis. 

ERROR  -  None. 

TUMBLE  rate 

USE  —  This  sets  up  the  parameter  value  for  modeling  inplane  tumbling  to  allow  com¬ 
putation  of  the  body  aspect  angle. 

RESPONSE  —  The  body  axis  is  established  in  the  trajectory  plane  and  is  given  a  tum¬ 
bling  rate  of  'rate'  dcg/sec. 

ERROR  —  The  body  orientation  may  be  specified  either  by  'POINTING'  or  by  'TUMBLE', 
but  not  by  both.  Whichever  is  more  recent  will  determine  which  is  used. 

SIGMA_TABLE  i,  number 

USE  —  This  command  is  the  means  by  which  a  cross-section  table  for  an  axisymmetric 
body  is  defined  to  allow  computation  of  cross  section  by  linear  interpolation. 

RESPONSE  —  This  command  puts  DEANE  in  a  mode  to  receive  a  list  of  'number' 

(Gj,  cr(G.))  pairs  which  define  cross-section  table  1  if  starting  on  the  next  input  line. 

These  pairs  should  be  entered  just  as  a  sequence  of  numbers  0^,  a(01),  G  ,  cr(G^),  .  .  . 
and  they  may  extend  over  several  input  lines,  a  (Q ^ )  is  interpreted  as  square  meters. 
This  tabic  naturally  overwrites  any  previous  cross-section  table  of  the  same  index  'i'. 

ERROR  —  There  is  an  upper  limit  of  20  pairs  and  'number'  must  be  an  integer.  The 
list  must  be  ordered  beginning  with  0°  and  terminating  with  (7(180°).  The  angles  must 
be  between  0°  and  180°,  and  ’i’  must  be  between  1  and  5. 

RV_TYPE  i 

USE  —  This  command  specifies  which  of  the  cross-section  tables  defined  is  to  be  used 
for  attaining  cross  sections. 

RESPONSE  —  The  cross-section  table  to  be  used  is  noted  as  cross-section  table  ' i' . 
ERROR  —  'i'  must  be  an  integer  between  1  and  5. 

RADAR_TYPE  n,  a,  b,  c,  0,  f,  E,  L,  T 

USE  —  For  determining  the  S/N  sensor,  parameters  need  to  be  known;  and  for  deter¬ 
mining  the  geometric  limits  as  to  what  the  sensor  can  see  (for  purposes  of  launching 
interceptors,  for  example),  additional  parameters  need  to  be  specified. 

'RADAR__TYPE'  achieves  this  by  associating  these  values  with  a  number  indicating  a 
radar  type,  and  the  user  can  then  associate  a  radar  type  with  each  radar  in  a  deployment. 

RESPONSE  —  Radar  type  'n'  is  defined  to  be  limited  by  'a'  kilometers  range,  ±'b'  de¬ 
grees  azimuth  off  azimuth  boresight,  ' c 1  degrees  lower  elevation  limit,  and  to  have  a 
beamwidth  of  ' Q '  degrees,  frequency  1  f '  megahertz,  'E*  joules  radiated  energy  per 
pulse,  'L'  decibel  loss,  and  an  operating  temperature  of  'T'  degrees  Kelvin.  (A  filled 
aperture  is  assumed  for  the  S/N  calculation.)  Type  1  is  reserved  for  a  default  radar 
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having  5500-km  range,  ±180°  azimuth,  1°  horizon,  1°  beamwidth,  500-MHz  frequency, 
125,000  joules/pulse,  4-dB  loss,  and  an  operating  temperature  of  1000 °K. 

ERROR  —  'n'  must  be  an  integer  between  2  and  5. 

NEW_SENSORS 

USE  —  This  command  allows  a  list  of  sensors  to  be  defined. 

RESPONSE  —  DEANE  enters  a  mode  to  receive  a  list  of  sensor  mnemonics,  latitudes, 
longitudes,  heights,  azimuth  boresights,  and  types.  Each  set  of  data  defining  a  sensor 
is  entered  on  one  line  (in  the  order  above)  and  the  list  is  terminated  by  typing  the  mne¬ 
monic  ’END1.  Default  value  for  the  height  is  0.0  km,  for  the  boresight  is  0.0°,  and  for 
the  type  is  1  if  the  numerical  entries  on  a  line  are  exhausted  before  satisfying  all  the 
values.  Any  previous  sensor  list  is  destroyed  at  the  beginning  of  this  command. 

ERROR  —  The  'type'  must  be  an  integer  between  1  and  5.  If  the  ’mnemonic'  entered 
for  the  addition  is  already  being  used  for  another  sensor,  an  error  message  will  say 
so  and  that  input  line  will  be  ignored.  Limits  of  12  characters  for  the  mnemonics  and 
a  total  of  40  sensors  arc  enforced. 

USE_SENSORS 

USE  —  This  command  selects  which  sensors  from  the  total  sensor  list  are  to  be  used 
either  for  purposes  of  observables  printout,  for  aiding  in  conducting  an  intercept,  or 
for  ” cookie-cutting"  functions.  It  also  allows  for  some  degree  of  editing  the  list  of 
sensors. 

RESPONSE  —  DEANE  enters  a  mode  to  receive  a  list  of  sensor  mnemonics,  one  per  line. 
This  list  is  terminated  by  entering  the  mnemonic  ’END1.  If  the  mnemonic  is  found  in 
the  existing  sensor  list,  that  sensor  is  marked  as  in  use.  If  a  previously  undefined 
mnemonic  is  followed  by  a  latitude,  longitude,  and,  optionally,  height,  boresight,  and 
type,  the  sensor  so  specified  will  be  added  to  the  existing  sensor  list  without  changing 
any  of  the  previous  entries,  and  this  sensor  will  be  marked  as  in  use.  If  a  previously 
defined  sensor  mnemonic  is  followed  by  a  latitude,  longitude,  and,  optionally,  height, 
boresight,  and  type,  the  said  values  for  the  named  sensor  will  be  changed  to  the  values 
entered  and  the  sensor  will  be  marked  as  in  use.  If  the  mnemonic  entered  is  ’ALL’, 
all  the  sensors  in  the  sensor  list  will  be  marked  as  in  use  and  the  command  will  be 
terminated. 

ERROR  —  There  is  a  limit  of  12  characters  for  each  mnemonic  and  a  total  of  40  sensors. 
Any  ’type’  entered  must  be  an  integer  between  1  and  5.  If  the  mnemonic  is  not  found 
in  the  sensor  list  and  is  not  followed  by  defining  data,  a  comment  to  that  effect  is  printed. 

TIMES  t,  t~  dt 
1  2 

USE  —  This  is  the  command  which  produces  the  printout  of  the  radar  observables. 

RESPONSE  —  For  each  sensor  listed  in  'USE_SENSORS',  values  for  the  quantities  listed 
in  'PRINT'  are  printed  using  the  working  trajectory.  Each  line  of  output  corresponds 
to  a  specific  time,  and  the  time  interval  covered  by  the  printout  is  from  1 1 ^ 1  to  ' t^'  in 
increments  of  'dt'.  The  calculation  of  the  body  aspect  angle  is  determined  by  either 
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1  -3-13053] 


launch  70  120  300 
Impact  39  -77  300 
parameter  1  4  304  1 
kepi e  r 

TOTAL  FLIGHT  TIME  IS  6475.12  SECONDS. 
ORBIT  ATTAINED  WITH  PERIOD  =  5406.176 

orint  time  lat  lonp  altitude  end 
times  0  -210  -50 

NO  SENSORS  ARE  IN  USE. 
u  se__senso  r  s 
Washington  39  -77 
end 

times  0  -210  -50 

OUTPUT  FOR  SENSOR  WASHINGTON 


TIME 

ALTITUDE 

LAT 

LONG 

0.0 

299.99 

39.00 

-77. 00 

-50.0 

300.68 

42.  30 

-76.51 

-100.0 

301.31 

45.61 

-75.98 

-150.0 

301.88 

48.91 

-75.42 

-200.0 

302.39 

52.21 

-74.81 

imes  6275.1  6576  50 

OUTPUT 

FOR  SENSOR 

WASHINGTON 

TIME 

ALTITUDE 

LAT 

LONG 

6275.1 

302.39 

52.21 

-74.81 

6325.1 

301.88 

48.91 

-75.42 

6375.1 

301.31 

45.61 

-75.98 

6425.1 

300.69 

42.31 

-76.51 

6475.1 

299.99 

39.00 

-77.00 

6525.1 

299.24 

35.70 

-77.47 

6575.1 

298.43 

32.39 

-77.91 

Fig.  I\  -1.  DEANE  commands  for  Example  4. 
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’TUMBLE’  or  ’POINTING',  depending  on  which  one  is  more  recent.  The  body  cross 
section  is  obtained  by  a  linear  interpolation  in  the  cross-section  tabic  of  type  ’RV  TYFT 
The  signal-to-noise  ratio  is  determined  using  the  parameters  describing  the  radar  type 
of  each  sensor  in  use  as  specified  via  1  RAL)AR_T YPE' .  For  this  command,  the  times 
may  be  positive  or  negative.  A  minus  time  is  time  to  impact,  and  a  positive  time  is 
time  after  launch  for  ’KEPLER’  and  'INTEGRATE'  and  time  after  ejection  for  'EJECT'. 

’ t ^ ’  and  ’t^'  may  stand  in  any  comparative  relationship. 

ERROR  —  This  command  is  terminated  if  either  of  the  following  situations  exists: 

(1)  Either  BAANGLE,  SIGMA,  or  S/N  is  requested  by  ’PRINT’ 
and  orientation  is  determined  by  'POINTING',  but  the  altitude 
set  in  'POINTING'  is  greater  than  apogee. 

(2)  Either  SIGMA  or  S/N  is  requested  by  'PRINT'  but  the  cross- 
section  table  type  specified  by  'R  V_T  YPE'  has  not  been  defined. 

The  output  for  any  particular  sensor  is  not  generated  if  the  parameters  for  that  radar 
type  have  not  yet  been  specified  via  'RADAR_TYPE'  and  S/N  is  requested  by  'PRINT'. 

Example  4 

To  see  how  the  positive  time  convention  is  used,  consider  the  set  of  DEANE  commands 
shown  in  Fig.  1V-1  where  the  first  four  commands  define  an  orbital  trajectory  which  passes 
over  the  impact  point  after  one  complete  period.  The  'PRINT'  command  states  what  is  to  be 
printed,  and  a  first  attempt  at  'TIMES'  yields  an  error  message.  Even  though  all  the  data  re¬ 
quested  is  sensor  independent,  some  sensor  must  be  defined  as  in  use.  This  is  done  via  the 
next  three  lines  where  the  editing  feature  of  ' USE__SENSORS'  is  used  to  define  and  mark  as  in 
use  a  sensor  at  the  impact  point.  The  execution  of  'TIMES  0  -210  -50’  gives  the  requested  out¬ 
put.  The  next  execution  of  'TIMES'  gives  equivalent  information  but  in  a  different  time  order 
and  with  the  positive  time  convention.  The  time  6475.1  is  equivalent  to  time  0,  the  time  6275.1 
is  equivalent  to  time  -200,  etc.  Note  that  data  beyond  the  "impact”  point  may  be  obtained  by 
using  the  positive  time  convention  as  seen  in  this  second  execution  of  'TIMES'.  This  use  of 
positive  time  is  the  motivation  for  printing  the  flight  time  for  'KEPLER'  and  'EJECT'. 

Example  5 

Probably  the  most  frequent  type  of  calculation  desired  is  simply  the  trajectory  observables 
for  a  specified  sensor  deployment  and  trajectory  geometry.  Consider  the  deployment  of  two 
acquisition  radars  ~1000km  apart  near  the  northern  border  of  the  continental  United  States  and 
an  RV  cross-section  model  similar  to  that  shown  in  Fig.  IV-2.  (The  radar  range  at  which  detec¬ 
tion  is  made  depends  on  which  cross  section  is  seen  by  the  radars.)  The  question  is:  "Can  such 
an  RV  on  a  minimum  energy  trajectory  from  China  targeted  at  one  radar  be  oriented  such  that 
detection  by  either  radar  will  be  made  on  the  smaller  cross  section  and,  hence,  later  in  time 
than  if  detection  could  be  made  on  the  larger  cross  section?"  Assume  that  the  radar  has  a  1  ° 
horizon  limit  and  the  sensitivity  and  detection  logic  are  such  that  detection  will  occur  at  2000-km 
range  on  the  small  cross  section. 

The  set  of  DEANE  commands  and  their  associated  output  seen  in  Fig.  IV- 3  indicate  that  such 
an  orientation  does  exist.  The  first  six  lines  in  the  figure  define  the  sensors  and  the  next  four 
lines  define  the  trajectory.  The  cross-section  table  is  then  defined  and  associated  with  the  RV 
via  'RV  TYPE’.  Interesting  observables  are  specified  in  'PRINT'  and  then  the  orientation  is 
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Fig.  IV- 3.  DEANE  commands  for  Example 


defined  in  'POINTING’.  Printout  from  'TIMES'  then  shows  that  first  detection  will  be  made  on 
the  small  cross  section  at  2000-km  range  by  the  radar  under  attack. 

ALTITUDE  alt 

USE  —  This  command  makes  available  to  the  user  the  time  at  which  the  object  on  the 
working  trajectory  is  at  a  certain  altitude.  It  is  intended  as  a  convenience  for  selecting 
times  of  interest  for  input  to  'TIMES’. 

RESPONSE  —  The  time  at  which  the  object  on  the  working  trajectory  is  at  altitude  'alt' 
is  printed.  If  'alt'  is  positive,  this  is  taken  at  the  impact  side  of  apogee;  if  'alt' 
is  negative,  this  is  taken  as  | ' alt*  |  kilometers  altitude  on  the  launch  side  of  apogee. 

ERROR  —  If  | 'alt'  |  is  greater  than  the  trajectory  apogee  or  if  |  ’ alt ’  |  is  lower  than  the 
trajectory's  minimum  altitude,  a  message  to  this  effect  is  printed  and  the  command 
is  terminated. 

VISIBLE  mnemonic 

USE  —  This  command  makes  available  to  the  user  the  earliest  and  latest  times  at  which 
a  sensor  can  view  a  trajectory.  Like  'ALTITUDE1,  this  is  intended  as  a  convenience 
for  selecting  input  times  for  ’TIMES'. 

RESPONSE  —  'mnemonic'  is  taken  as  the  name  of  a  sensor  in  the  sensor  list  and  it 
accordingly  has  a  position,  boresight,  and  a  limiting  azimuth,  elevation,  and  range 
as  specified  through  the  sensor  type.  This  is  sufficient  information  to  calculate  the 
times  at  which  this  sensor  can  first  and  last  view  an  object  on  the  working  trajectory. 

If  the  sensor  can  view  the  trajectory,  these  times  and  the  associated  ranges,  azimuths, 
and  elevations  are  printed.  If  the  sensor  cannot  view  the  trajectory,  a  comment  is 
printed  to  indicate  this.  If  'mnemonic'  is  'ALL’,  the  above  information  is  printed  for 
all  the  sensors  marked  as  in  use. 

ERROR  —  If  ’mnemonic'  is  not  found  in  the  sensor  list  and  is  not  'ALL',  a  message 
will  so  indicate  and  the  command  will  be  terminated.  Also,  if  the  corresponding  sensor 
parameters  are  undefined,  a  message  will  so  indicate. 

Example  6 

To  see  how  'VISIBLE'  can  be  used,  in  Example  5  we  needed  observables  printout  over  some 
time  interval.  With  the  trajectory  and  radar  geometry  of  Example  5,  ’VISIBLE’  may  be  used 
as  seen  in  Fig.  IV-4  to  indicate  that  the  RV  first  crosses  the  radar  horizon  at  —592  sec.  Thus, 

—  600  is  a  good  time  to  start  requesting  observables.  If  the  radar  had  a  2°  horizon  instead  of  the 
1°  default  horizon,  we  would  have  to  define  this  via  'RADAR  TYPE'  as  seen  in  Fig.  IV-4.  . 
'VISIBLE'  then  indicates  a  horizon  crossing  at  —560  sec. 
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V.  TRAJECTORY  INTERCEPTION  COMMANDS 


FLYOUT_TABLE  i,  dela,  delr,  numa,  numr 

USE  -  This  command  provides  the  means  for  defining  interceptor  flyout  characteristics. 

RESPONSE  —  The  flyout  reach  table  to  be  entered  should  resemble  that  shown  in  Fig.  V-l, 
where  the  entries  are  the  flyout  times  in  seconds  required  to  make  an  intercept  at  the 
specified  tangent  range  from  the  farm  and  at  the  specified  altitude  above  the  earth.  For 
(tangent  range,  altitude)  pairs  the  interceptor  cannot  reach,  1.0E7  should  be  entered  for 
the  flyout  time.  The  interceptor  type  is  'i',  'dela'  is  the  table  altitude  increment  in 
kilometers,  'delr'  is  the  table  tangent  range  increment  in  kilometers,  'numa'  is  the 
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ALTITUDE 

TANGE  NTX. 
RANGE  \ 

0 

dela 

2  x  dela 

(numo-1)  x  dela 

0 

0,0 

*0,1 

*0,2 

*0,  numa-1 

delr 

Vo 

t 

Li 

*1,2 

*1,numo  -1 

2  X  delr 

( numr  -1)  x  delr 

*numr-1,0 

*numr-1,  numa-1 

Fig.  V-l.  General  form  of  interceptor  reach  table. 

number  of  altitude  entries,  and  'numr'  is  the  number  of  tangent  range  entries.  After 
reading  these  numbers,  DEANE  expects  a  sequence  of  'numa'  x  'numr'  numbers  starting 
on  the  next  input  line  to  correspond  to  the  table  entries.  The  order  of  the  entries  is 
row  by  row,  i.e.,  the  sequence  would  be  tn  n,  tn  tn  _  , ,  t.  n,  .  .  .  .  These* 

quence  may  extend  over  any  number  of  input  lines,  and  it  overwrites  any  previous  flyout 
table  of  type  'if. 

ERROR  -  'i'  must  be  an  integer  between  1  and  5,  and  'numa'  and  'numr'  must  be  integers 
no  greater  than  30. 

NEW  FARMS 

USE  -  This  command  allows  a  list  of  interceptor  farms  to  be  defined. 

RESPONSE  -  DEANE  enters  a  mode  to  receive  a  list  of  farm  mnemonics,  latitudes, 
longitudes,  heights,  populations,  and  types.  Each  set  of  data  defining  a  farm  is  entered 
on  one  line  in  the  order  designated  above,  and  the  list  is  terminated  by  typing  the  mne¬ 
monic  'END'.  Default  value  for  the  height  is  0km,  for  the  population  (i.e.,  number  of 
interceptors  at  the  farm)  is  0,  and  for  the  type  is  1  if  the  numerical  entries  of  a  line 
are  exhausted  before  satisfying  all  the  values.  Any  previous  farm  list  is  destroyed  at 
the  beginning  of  this  command. 

ERROR -The  'type'  must  be  an  integer  between  1  and  5.  If  the  'mnemonic'  entered  for 
the  addition  is  already  being  used  for  another  farm,  an  error  message  will  say  so  and 
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that  input  line  will  be  ignored.  Limits  of  12  characters  for  mnemonics  and  a  total  of 
40  farms  are  enforced. 

USE_FARMS 

USE  -  This  command  selects  which  farms  from  the  total  farm  list  are  to  be  used  as  the 
origin  of  interceptors.  It  also  permits  editing  the  farm  list  to  some  extent. 

RESPONSE  -  DEANE  enters  a  mode  to  receive  a  list  of  farm  mnemonics,  one  per  line, 
and  this  list  is  terminated  by  entering  the  mnemonic  'END'.  If  the  mnemonic  is  found 
in  the  farm  list,  that  farm  is  marked  as  in  use.  Should  an  undefined  mnemonic  be  fol¬ 
lowed  by  a  latitude,  longitude,  and,  optionally,  height,  population,  and  type,  the  farm 
so  specified  will  be  added  to  the  farm  list  without  changing  any  of  the  previous  entries 
and  it  will  be  marked  as  in  use.  If  a  defined  farm  is  followed  by  a  latitude,  longitude, 
and,  optionally,  height,  population,  and  type,  the  said  values  for  that  farm  will  be 
changed  to  the  values  entered  and  the  farm  will  be  marked  as  in  use.  If  the  mnemonic 
entered  is  'ALL',  all  the  farms  in  the  farm  list  will  be  marked  as  in  use  and  the  com¬ 
mand  will  be  terminated. 

ERROR  —  There  is  a  limit  of  12  characters  for  each  mnemonic  and  a  total  of  40  farms. 
Any  ’type'  entered  must  be  an  integer  between  1  and  5.  If  the  mnemonic  is  not  found  in 
the  farm  list  and  is  not  followed  by  defining  data,  a  comment  to  that  effect  is  printed. 

T  R  A  CK1NG_SE  NSOR  S 

USE  -  This  command  defines  which  sensors  are  to  be  used  for  modeling  tracking  to 
intercept. 

RESPONSE  -  DEANE  enters  a  mode  to  receive  a  list  of  sensor  mnemonics,  one  per  line, 
terminated  by  the  mnemonic  'END'.  Each  sensor  entered  is  then  noted  as  a  site  from 
which  tracking  can  be  done. 

ERROR  -  If  any  mnemonic  entered  is  not  in  the  sensor  list,  a  message  to  that  effect  is 
printed  and  that  mnemonic  is  ignored. 

RESTRICT  low,  high,  accuracy,  timel,  time2 

USE  -  This  command  specifies  values  for  restrictions  to  be  placed  on  intercepts. 

RESPONSE  -  ’low'  is  the  value  of  the  minimum  intercept  altitude,  'high'  is  the  altitude 
below  which  the  object  must  be  when  the  interceptor  is  launched  (i.e.,  maximum  com¬ 
mit  altitude),  and  'accuracy'  is  the  time  accuracy  desired  in  determining  the  interceptor 
launch  times,  'timel'  is  the  required  time  delay  before  the  interceptor  launch  after 
the  first  possible  detection  by  any  sensor  marked  as  in  use,  and  'time2'  is  the  length  of 
the  time  interval  immediately  prior  to  the  intercept  during  which  the  object  must  be 
visible  to  some  sensor  being  used  for  tracking. 

ERROR  —  All  entries  must  be  positive,  and  'low'  must  be  less  than  'high'. 

INTERCEPT 

USE  -  This  command  provides  information  concerning  intercepts  of  an  object  on  the 
working  trajectory. 
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RESPONSE  -  The  times  of  the  intercepts  with  the  earliest  intercept  time  and  with  the 
latest  interceptor  launch  time  are  printed  together  with  the  farms  used  and  the  intercept 
altitudes  and  ranges  from  the  farms.  The  data  defining  acceptable  intercepts  are  spec¬ 
ified  via  ’  RESTRICT',  ' USE_SENSORS',  ' RADAR_T YPE',  'USE_FARMS\  ’FLYOUT_ 
TABLE',  and  ’  TRACK1NGJ3ENSORS'.  ' USE_FARMS'  defines  a  set  of  farms  from  which 
interceptors  can  be  launched;  the  type  of  interceptor  at  a  farm  references  a  set  of  fly¬ 
out  characteristics  specified  via  '  FLYOUT_TABLE',  which  then  define  the  performance 
of  interceptors  launched  from  that  farm.  Intercepts  must  occur  at  or  before  the  min¬ 
imum  intercept  altitude  specified  in  'RESTRICT',  and  interceptor  launches  must  occur 
at  or  after  the  maximum  commit  altitude  set  in  'RESTRICT'.  If  there  are  no  sensors 
marked  as  in  use  and  there  are  no  tracking  sensors  defined,  the  above  constitute  the 
only  restrictions  on  intercepts.  If  there  are  some  sensors  in  use,  before  any  intercep¬ 
tor  launch  occurs  some  sensor  of  those  marked  in  use  via  'USE_SENSORS'  must  have 
been  able  to  detect  the  object  as  limited  by  the  elevation,  azimuth,  and  detection  range 
associated  with  those  sensors  via  'RADAR_TYPE'.  The  earliest  acceptable  interceptor 
launch  is  then  'timel'  sec  after  this  radar  detection  where  'timel'  is  set  in  'RESTRICT'. 
This  allows  modeling  the  time  needed  to  insure  acquisition  in  the  surveillance  pattern 
plus  the  time  needed  to  establish  a  track  sufficiently  good  to  permit  committing  an  in¬ 
terceptor.  If  there  are  some  sensors  being  used  for  tracking  as  specified  in  'TRACK- 
ING_SENSORS\  one  of  these  sensors  must  be  able  to  track  to  intercept  as  modeled  by 
being  able  to  view  the  object  during  the  'time2'  sec  just  prior  to  the  intercept,  where 
'time2f  is  set  in  'RESTRICT1. 

If  no  acceptable  intercepts  can  be  made,  a  4-bit  code  is  printed  to  give  insight  into  the 
situation.  The  first  (leftmost)  bit  is  0  if  the  earliest  launch  is  determined  by  the  spec¬ 
ified  maximum  commit  altitude,  and  it  is  1  if  the  acquisition/track  restriction  limits 
the  earliest  launch.  The  next  bit  is  0  if  the  lowest  intercept  altitude  is  determined  by 
the  minimum  intercept  altitude  set  in  'RESTRICT',  and  is  1  if  it  is  determined  by  the 
tracking  to  intercept  restriction.  The  third  bit  is  0  if  an  interceptor  cannot  physically 
reach  the  trajectory  between  the  two  limits  set  above.  If  one  interceptor  can  reach  the 
trajectory  but  the  flight  time  is  too  long  to  be  an  acceptable  intercept,  the  third  bit  is  1. 
If  the  final  bit  is  0,  this  means  that  no  intercept  can  be  made  with  intercept  above  the 
lowest  acceptable  altitude  and  launch  after  the  earliest  acceptable  time;  if  the  final  bit 
is  1,  such  intercepts  can  be  made  but  are  not  acceptable  due  to  a  lack  of  the  required 
visibility  modeling  track  to  intercept.  For  example,  the  code  "0000"  would  indicate 
that  the  interceptor  does  not  have  enough  performance  to  make  any  intercepts  within  the 
launch-intercept  altitude  band  specified  in  'RESTRICT';  the  code  "4  010"  would  indicate 
that  the  minimum  intercept  altitude  set  in  'RESTRICT'  specifies  the  lowest  acceptable 
intercept  altitude  and  intercepts  can  be  made  above  this  altitude,  but  none  can  be  made 
which  have  a  launch  after  the  earliest  acceptable  launch  time  as  determined  by  the 
acquisition/track  restriction. 

ERROR  -  If  no  farms  are  in  use  or  the  flyout  characteristics  for  any  farm  in  use  are 
not  defined,  a  message  will  so  indicate  and  the  command  will  be  terminated.  If  radar 
restrictions  are  included  and  the  radar  parameters  for  any  such  radar  are  not  defined, 
a  message  will  be  printed  to  this  effect  and  the  command  will  be  terminated. 
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Example  7 


Consider  the  case  of  a  defense  site  in  the  central  United  States  which  consists  of  a  radar 
and  an  interceptor  farm.  The  radar  must  provide  its  own  surveillance  and  the  threat  of  interest 
is  from  China  with  an  assumed  RV  cross  section  such  that  the  radar  ean  detect  it  at  4000-km 
range.  The  radar  has  a  1.5°  horizon  and,  because  of  a  narrow  threat  corridor,  we  assume  it 
has  full  azimuth  coverage.  (This  assumption  will  speed  execution.)  It  is  further  assumed  that 
60-sec  time  delay  after  first  possible  geometric  detection  will  give  guaranteed  detection  and 
enough  initial  tracking  to  permit  an  interceptor  launch.  The  (idealized)  exoatmospheric  intercep¬ 
tor  at  the  site  has  the  reach  curves  shown  in  Fig.  V-2,  and  it  has  a  minimum  intercept  altitude 
of  7  5  km  due  to  warhead  considerations.  What  we  want  to  know  is,  assuming  threats  on  the  radar, 
when  nan  intercepts  be  made? 


Fig.  V-2.  Interceptor  reach  curves 
for  Example  7. 


z  4  e  s  10  iz 
TANGENT  RANGE  (km  x  10  2 ) 

The  DEANE  commands  seen  in  Fig.  V-3  provide  the  answers  to  this  question.  (The  is 
used  as  a  break  character  to  mark  the  start  of  a  new  DEANE  line  for  processing  while  staying  on 
the  same  teletype  line.)  The  first  line  in  Fig.  V-3  defines  the  radar  deployment  with  pertinent 
characteristics,  and  the  second  line  defines  the  farm  deployment.  The  next  five  lines  define  the 
interceptor  performance,  and  then  it  is  specified  that  there  are  no  tracking  to  intercept  radars. 
(In  this  geometry,  the  tracking  to  intercept  visibility  requirement  would  not  be  a  limiting  factor 
so  we  neglect  it  to  speed  execution.)  Now  the  restrictions  are  set.  The  minimum  intercept  alti¬ 
tude  is  set  to  7  5  km  and  the  maximum  commit  altitude  is  set  very  large  so  that  the  radar  detec¬ 
tion  will  be  limiting.  A  2- sec  accuracy  is  sufficient,  and  the  time  delay  of  60  sec  is  entered. 

With  all  this  necessary  data  defined,  the  next  line  defines  and  calculates  a  minimum  energy 
trajectory  impacting  at  the  site.  'INTERCEPT1  causes  the  next  four  lines  of  output.  Because 
of  the  geometry  here,  the  latest  launch  gives  the  latest  intercept  which  is  at  7  5-km  altitude. 
Execution  of  'VISIBLE'  shows  that  the  earliest  intercept  corresponds  to  interceptor  launch  60  sec 
after  first  geometric  visibility.  It  is  interesting  to  note  that  the  latest  launch  occurs  before  the 
earliest  intercept  so  that  the  system  could  not  attempt  one  intercept,  evaluate  the  intercept,  and 
then  launch  a  second  interceptor  ( i. e. ,  this  system  has  no  shoot- look- shoot  capability  against  a 
minimum  energy  trajectory).  The  trajectory  is  then  changed  to  10°  re-entry  and  'INTERCEPT' 
says  that  no  intercepts  can  be  made  and  the  code  is  "1010"  (discussed  in  command  definition  for 
'INTERCEPT').  This  is  understandable  for  two  reasons:  first,  detection  will  be  horizon  limited 
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200 

parameter  0 

2  20 

contour 

10  180  25  - 

20 

0 

THE  CENTER 

IS  37 

.000  LAT  AND 

-97. 000 

LONG. 

PSI 

RANGE 

LAT 

LONG 

CODE 

0.0 

425.0 

40.82 

-97. 00 

1010 

10.0 

425.0 

40.76 

-96.12 

1010 

20.0 

450.0 

40.79 

“95.17 

1010 

30.0 

450.0 

40.47 

“94.34 

1010 

40.0 

475.0 

40.  22 

“93.41 

1010 

50.0 

525.0 

39.94 

“92.29 

1010 

60.0 

575.0 

39.44 

“91. 20 

1010 

70.0 

625.0 

38.73 

“90.23 

1010 

80.0 

725.0 

37.86 

“88.86 

1010 

90.0 

825.0 

36.64 

“87. 75 

1010 

100.0 

950.0 

35.07 

-86.71 

1010 

110.0 

1025.0 

33.39 

“86.62 

0000 

120.0 

1100.0 

31.62 

-86.94 

0000 

130.0 

1200.0 

29.68 

-87.50 

0000 

140.0 

1375.0 

27.19 

-88.10 

0000 

150.0 

1700.0 

23.  47 

-88. 74 

0000 

160.0 

2150.0 

18.64 

-90.14 

0000 

170.0 

2625.0 

13.69 

-92.90 

0000 

180.0 

3125.0 

8.91 

-97.00 

0000 

Fig.  V-4.  DEANE  commands  for  Example  8. 
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and  will  occur  later  in  time  than  for  the  minimum  energy  case;  second,  75-km  altitude  occurs 
farther  uprange  for  a  10°  trajectory  than  for  a  minimum  energy  trajectory.  We  thus  have  farther 
to  go  in  less  time  than  in  the  minimum  energy  case.  The  trajectory  is  now  changed  to  30°  re-entry 
and  we  note  that,  we  do  have  a  shoot-look- shoot  capability  with  the  earliest  intercept  at  -122  sec 
and  the  latest  launch  at  -108  sec. 

CONTOUR  dclpsi,  psi,  dr,  tau,  alt,  beta 

USE  -  This  command  produces  a  table  defining  footprints  of  interceptor  coverage. 

RESPONSE  -  A  contour  of  ground  coverage  is  printed  in  the  form  of  an  angle  and  ground 
range  from  some  center  point,  which  is  also  printed.  The  corresponding  latitudes  and 
longitudes  are  printed  as  is  a  code  indicating  the  limiting  factor  for  that  point  of  the  con¬ 
tour.  The  limitations  are  the  same  as  those  discussed  for  the  command  'INTERCEPT' 
and  have  the  same  meaning.  The  interpretation  of  the  contour  is  that  if  the  impact  point 
is  'alt'  km  above  a  point  interior  to  the  contour  (using  the  launch  site  set  in  '  LAUNCH' 
and  trajectory  parameters  set  in  'PARAMETER'),  the  threatening  object  can  be  inter¬ 
cepted  using  the  same  intercept  restrictions  discussed  in  'INTERCEPT'.  Similarly,  a 
point  exterior  to  the  contour  cannot  be  defended  under  the  same  restrictions.  If  'beta' 
is  specified  and  is  nonzero,  the  trajectories  used  to  obtain  the  contour  are  Keplerian 
with  the  re-entry  portion  obtained  by  integration  using  a  real  atmosphere  and  the  ballis¬ 
tic  coefficient  history  specified  in  '  BETA  _TABLE!.  If  the  contour  is  desired  for  strictly 
Keplerian  trajectories,  'beta'  should  be  omitted  from  the  input. 

The  table  generated  is  in  angular  increments  of  'delpsi'  deg  from  zero  to  'psi'  deg  with 
a  radial  accuracy  of  'dr'  km.  If  'tau'  is  less  than  zero,  the  contour  is  for  single  inter¬ 
cepts;  but  if  'tau'  is  greater  than  zero,  the  contour  is  for  shoot-look- shoot  with  'tau'  see 
between  the  first  intercept  and  the  second  launch. 

ERROR  -  The  errors  are  the  same  as  discussed  in  'INTERCEPT'.  If  the  contour's 
north-south  dimension  is  less  than  'dr',  a  message  will  so  indicate  and  the  command 
will  be  terminated.  A  word  of  warning  is  in  order  regarding  the  interpretation  of  the 
results  of  this  command  due  to  the  algorithm  used.  The  algorithm  finds  one  point  of  the 
curve  and  then  traces  the  rest  assuming  the  contour  is  a  single  valued  function  in  a  polar 
coordinate  system  with  the  origin  at  the  center  point  determined.  No  test  is  made  to 
find  undefended  areas  interior  to  the  contour. 

If  'beta'  is  entered,  it  must  be  an  integer. 

Example  8 

Consider  an  extension  of  the  problem  posed  in  Example  7.  With  the  same  defense  site  and 
characteristics,  we  want  to  know  what  the  defended  footprint  is.  Use  of  'CONTOUR'  in  the  se¬ 
quence  of  DEANE  commands  seen  in  Fig.  V-4  provides  this  information  for  a  20°  re-entry  tra¬ 
jectory  neglecting  tracking  to  intercept  considerations.  The  first  18  lines  in  Fig.  V-4  are  identical 
with  the  first  9  lines  in  Fig.  V-3,  but  the  launch  point  has  been  changed  and  a  nonrotating  earth 
is  used.  This  is  done  so  that  the  resulting  footprint  will  be  symmetric  about  the  farm  longitude, 
and  only  half  the  footprint  need  be  calculated.  It  is  interesting  to  interpret  the  codes  printed  in 
Fig.  V-4.  To  the  side  and  uprange  of  the  site,  the  footprint  is  limited  by  the  radar  acquisition 
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(CODE  1010).  To  the  rear,  the  footprint  is  limited  only  by  the  reach  capabilities  of  the 
interceptor  (CODE  -  0000).  The  results  of  'CONTOUR'  seen  in  Fig.  V-4  are  plotted  in  Fig.  V-5 
as  are  similar  results  for  10°  and  30°  re-entry  trajectories. 


THREAT 

\7 


1 18-3-13061 


Fig.  V-5.  Footprints  for  Example  8  and 
for  10°  and  30°  re-entry  trajectories. 


Example  q 

To  demonstrate  the  endoatmospheric  footprint  capability  with  an  integrated  trajectory,  con¬ 
sider  an  (idealized)  interceptor  which  can  fly  at  2  km/sec  with  a  maximum  time  of  flight  of 
60  sec.  Sited  with  a  radar,  we  assume  that  the  radar's  performance  is  described  by  a  maximum 
commit  altitude  and  that  the  radar's  field  of  view  does  not  geometrically  limit  intercepts.  The 
threat  of  interest  here  is  at  minimum  energy  with  a  0  of  1500  down  to  30  km,  and  a  /3  of  500 
from  30-  to  0-km  altitude.  With  a  minimum  intercept  altitude  of  6km  and  a  maximum  commit 
altitude  of  70  km,  the  DEANE  commands  in  Fig.  V-6  yield  the  contour  which  defines  the  defended 
region.  By  inspecting  the  codes,  note  that  the  contour  is  not  limited  by  the  interceptor  reach 
but  by  the  commit  and  minimum  intercept  altitudes.  This  footprint  is  plotted  in  Fig.  V-7,  with 
footprints  generated  for  50-  and  30-km  commit  altitudes.  Also  plotted  for  interest  is  the  foot¬ 
print  for  a  maximum  commit  altitude  of  30  km  and  a  minimum  intercept  altitude  of  3  km.  This 
footprint  illustrates  the  fact  that  there  is  significantly  more  slowdown  at  lower  altitudes  than  at 
the  higher  altitudes,  thus  providing  almost  as  much  coverage  as  for  commit  after  70  km  and 
intercept  by  6  km. 
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7^3- 13062  1 


f  1  you t__tab  1  e  5  20  20  7  7 

0  10  20  30  40  50  60 

10  14.1  22.4  31.6  41.2  51  1.0e7 

20  22.4  28.3  36  44.7  54  1 . Oe 7 

30  31.6  36  42.5  50  58  1.0e7 

40  41.2  44.7  50  57  1.0e7  1.0e7 

50  51  54  58  1 . 0e7  1.0e7  1.0e7 

60  1 . 0e7  1 . Oe  7  1.0e7  1.0e7  1.0e7  1.0e7 

use_ f a  rms 

farml  37  -97  0  0  5 
end 

use_sensors 

end 

track! ng_sensors 
end 

beta_table  1500  30  500 
parameter  0  1 
launch  38  83  200 
rest rict  6  70  0.5  0  0 
contour  10  180  1  -20  0  1 


THE  CENTER 

IS  37. 

000  LAT  AND 

-97. 000 

LONG. 

PSI 

RANGE 

LAT 

LONG 

CODE 

0.0 

76.0 

37.68 

-97.00 

0010 

10.0 

77.0 

37.68 

-96.85 

0010 

20.0 

77.0 

37.65 

-96.70 

0010 

30.0 

77.0 

37.60 

-96. 56 

0010 

40.0 

78.0 

37.54 

-96.43 

0010 

50.0 

79.0 

37.45 

-96.31 

0010 

60.0 

80.0 

37.36 

-96.22 

0010 

70.0 

81.0 

37.  25 

-96.14 

0010 

80.0 

82.0 

37.12 

-96.09 

0010 

90.0 

83.0 

37.00 

-96.07 

0010 

100.0 

84.0 

36.87 

-96.07 

0010 

110.0 

85.0 

36.  74 

-96.10 

0010 

120.0 

86.0 

36.61 

-96.17 

0010 

130.0 

87.0 

36.  50 

-96.25 

0010 

140.0 

88.0 

36. 39 

-96.37 

0010 

150.0 

89.0 

36.31 

-96.50 

0010 

160.0 

89.0 

36.25 

-96.66 

0010 

170.0 

90.0 

36.20 

-96.83 

0010 

180.0 

90.0 

36.19 

-97.00 

0010 

Fig.  V-6.  DEANE  commands  for  Example  9. 
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MINIMUM  ENERGY 
THREAT 


V7 


1  19-3-13063  1 


Fig.  V-7.  Footprint  for  Example  9  and  three 
other  endoatmospheric  cases. 
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VI.  "COOKIE -CUTTING"  COMMANDS 


SITES_ABOUT_SMSA  <p  d  rankl  rank2 

USE  -  This  command  allows  for  seeing  what  deployment  sites  (farms  or  radars)  are 
within  a  certain  distance  of  certain  SMSA's  (Standard  Metropolitan  Statistical  Areas). 

RESPONSE  -  For  SMSA’ s  of  'rankl*  through  'rank2',  the  rank,  name,  and  population 
of  the  SMSA  are  printed  under  which  are  printed  the  name  of  each  site,  the  azimuth  of 
the  SMSA  from  the  site,  and  its  distance  to  the  SMSA,  if  that  distance  is  not  greater 
than  'd'.  If  V  is  'FARMS',  the  farm  deployment  specified  by  'IJSE_F  ARMS'  is  used; 
if  V  is  'RADARS',  the  sensor  deployment  specified  by  'USESENSORS'  is  used  for  the 
calculation. 

ERROR  -  If  V  is  'FARMS'  and  no  farms  are  in  use  or  if  V  is  'RADARS'  and  no  radars 
are  in  use,  a  message  will  so  indicate  and  the  command  will  be  terminated,  'rankl' 
must  be  less  than  or  equal  to  'rank2',  and  the  range  of  ranks  is  from  1  to  228.  t<pi  must 
be  either  'FARMS'  or  'RADARS'. 

SMSAS_ABOUT_SIT  E  <p  name  d 

USE  -  This  allows  for  seeing  what  SMSA's  are  within  a  certain  distance  of  one  particular 
site. 

RESPONSE  -  For  each  SMSA  within  distance  ' d'  of  the  site  'name',  the  SMSA  rank, 
name,  population,  distance  to  the  site,  and  azimuth  from  the  site  are  printed.  If  V 
is  'FARM',  'name'  references  a  farm;  if  '<p*  is  'RADAR',  'name'  references  a  sensor. 

ERROR  -  If  the  site  'name'  is  not  in  the  list  of  farms  or  sensors  (as  indicated  by  V), 
a  message  will  so  indicate  and  the  command  will  be  terminated.  '<p'  must  be  either 
'FARM'  or  'RADAR'. 

IN  SID  E_SIT  E_RING  <p  dd  d2  rankl  rank2 

USE  -  This  command  considers  an  entire  farm  or  radar  deployment  and  determines 
which  SMSA's  are  in  a  certain  range  of  distances  from  any  site. 

RESPONSE  -  For  SMSA's  of  'rankl'  through  'rank2',  if  the  distance  from  the  SMSA  to 
the  closest  site  in  use  is  greater  than  ' d^ '  but  not  greater  than  '  d^'  the  SMSA  rank,  name, 
and  population  will  be  printed  with  the  name  of  the  closest  site,  that  distance,  and  the 
azimuth  of  the  SMSA  from  the  site.  If  V  is  'FARMS',  the  farm  deployment  specified 
by  'USE_FARMS'  is  used;  if  V  is  'RADARS',  the  sensor  deployment  specified  by 
'USEJSENSORS'  is  used  for  the  calculation. 

ERROR  -  If  no  sites  are  marked  in  use  (depending  on  V),  a  message  will  so  indicate 

and  the  command  will  be  terminated,  'd  '  must  be  less  than  'd  '  'rankl'  must  be  less 

1  2 

than  or  equal  to  'rank2',  and  the  range  of  ranks  is  1  to  228.  V  must  be  either  'FARMS' 
or  'RADARS'. 

CUM_DISTRIBUTION  <p  rankl  rank2 

USE  -  This  command  considers  a  deployment  of  farms  or  radars  and  prints  the  nor¬ 
malized  cumulative  population  distribution  as  a  function  of  distance  from  the  closest  site. 
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RESPONSE  -  The  rank,  name,  and  population  of  SMSA' s  'rankl'  through  'rank2'  are 
printed  in  order  of  increasing  distance  from  the  closest  site  in  use.  Also  printed  for 
each  SMSA  is  the  name  of  the  closest  site,  the  distance  to  that  site,  the  total  popula¬ 
tion  within  that  distance,  and  this  population  normalized  to  the  total  population  in  the 
SMSA's  of  'rankl'  through  'rank2'.  If  '<p'  is  'FARMS',  the  farm  deployment  specified 
by  'USE_FARMS'  is  used;  if  V  is  'RADARS',  the  sensor  deployment  specified  by 
'  USE_SENSORS'  is  used  for  the  calculation. 

ERROR  -  If  no  sites  are  in  use  (depending  on  1  <p'),  a  message  will  so  indicate  and  the 
command  will  be  terminated,  'rankl'  must  be  less  than  or  equal  to  'rank2',  with  the 
range  of  ranks  from  1  to  228.  V  must  be  either  'FARMS'  or  'RADARS'. 

INSIDE  BOTH  dd  d,  rank!  rank2 

USE  -  This  command  considers  a  deployment  of  both  farms  and  radars  and  determines 
what  SMSA's  are  within  a  certain  distance  of  a  farm  and  within  a  certain  distance  of 
a  radar. 

RESPONSE  —  For  SMSA's  of  'rankl'  through  'rank2',  if  that  SMSA  is  within  a  distance 
'd^'  of  any  farm  in  use  and  is  also  within  a  distance  'd^'  of  any  radar  marked  in  use, 
that  SMSA  rank,  name,  and  population  are  printed. 

ERROR  -  If  no  farms  are  in  use  or  no  radars  are  in  use,  a  message  will  so  indicate 
and  the  command  will  terminate,  'rankl'  must  be  less  than  or  equal  to  'rank2',  with 
the  range  of  ranks  being  from  1  to  228. 

NOT  INSIDE  BOTH  d,  d~  rankl  rank2 

USE  -  This  command  is  the  logical  complement  of  'INSIDE  BOTH'. 

RESPONSE  -  For  SMSA's  of  'rankl'  through  'rank2',  if  the  distance  from  the  SMSA 
to  the  nearest  farm  in  use  is  greater  than  ' d^ '  or  if  the  distance  to  the  nearest  radar 
in  use  is  greater  than  'd^',  the  SMSA  rank,  name,  and  population  are  printed.  Under 
this  are  printed  the  names  and  distances  to  the  closest  farm  and  the  closest  radar  if 
these  distances  are  greater  than  'd^'  and  1  d^ ' ,  respectively.  Note  that  by  setting  'd^' 
and  'd^'  to  zero,  the  result  is  a  listing  of  SMSA's  followed  by  the  name  and  distance 
to  the  closest  farm  in  use  and  the  name  and  distance  to  the  closest  radar  in  use. 

ERROR  -  If  no  farms  are  in  use  or  no  radars  are  in  use,  a  message  will  so  indicate 
and  the  command  will  be  terminated,  'rank!'  must  not  be  greater  than  'rank2',  with 
the  range  of  ranks  being  from  1  to  228. 
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Fig.  VI-1.  Population  defense  radar  deployment  for  Example  10. 


Example  10 

Consider  the  12- site  population  defense  radar  deployment  shown  in  Fig.  VI- 1.  The 
DEANE  commands  seen  in  Fig.  VI-2  provide  data  for  making  comments  about  this  deployment. 
'USE_SENSORS'  in  Fig.  VI-2  defines  the  radars  preparatory  to  the  cookie-cutting  commands. 

If  a  radar’s  coverage  extends  to  360  km  (200 nmi)  from  the  site,  then  ’SITES_ABOUT_SMSA' 
shows  that  three  of  the  top  five  SMSA’s  have  double  coverage.  (For  example.  New  York  City 
is  at  a  distance  of  172km  and  an  azimuth  of —158°  from  RADAR2,  and  at  a  distance  of  268km 
and  an  azimuth  of  7  3°  from  RADAR8. )  To  see  what  SMSA's  are  capable  of  defense  from 
RADAR12,  we  execute  ’SMSAS_ABOUT_SITE\  The  coverage  for  the  top  ten  SMSA's  is  found 
by  executing  'INSIDEJ3IT EARING',  and  we  note  that  all  the  top  ten  SMSA's  are  within  150  nmi 
(27  0  km)  from  a  site  and  none  are  closer  than  65  nmi  to  a  SMSA.  'CUM_DISTRIBUTION'  gives 
coverage  information  for  the  next  fifteen  largest  SMSA's,  and  we  see  that  all  but  two  are  within 
150  nmi  of  a  site  and  all  are  within  200  nmi  of  a  site. 
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TT3-<3065| 


usp.sensors 

radarl  47  -122  #  radar2  42.1  -73.2  #  radar3  33.8  -115.4 
radar4  31.5  -97.3  *  radar5  41.3  -84.7  #  radar6  32.6  -87.8 
radar7  28.1  -80.9  #  radar8  40  -77  #  radar9  38  -121 
radarlO  44.1  -90.8  #  radarll  38.5  -93  #  radarl2  40  -108 


end 

s i tes_ 

about  smsa  radars 

360  1  5 

1 

NEW  YORK, NY 

10694632. 

RADAR2 

171.79  -157.9 

RADAR8 

267.72  72.9 

2 

CHICAGO, 1 LL 

6220913. 

RADAR5 

260.72  -75.8 

RADAR  10 

353.92  134.4 

3 

LOS  ANGELES, CAL 

6038771. 

RADAR3 

259.56  -84.3 

4 

PHI LADELPHI A, PA 

4342897. 

RADAR2 

296.04  -146.4 

RADAR8 

161.32  94.7 

5 

DETROIT, HIGH. 

3762360. 

RADAR  5 

164.86  47.1 

smsas 

about  s i te  rada r 

radarl2  360 

RANK 

NAME 

POPULATION 

DISTANCE  AZIMUTH 

27 

DENVER, COLO 

929383. 

257.  66 

95.2 

63 

SALT  LAKE  CITY, UTAH  447795. 

342. 32 

-74.6 

164 

COLORADO  S PROS , COLO  143742. 

301. 60 

114  .  5 

188 

PUEBLO, COLO 

118707. 

347. 71 

122.  3 

196 

PROVO, UTAH 

106991. 

313.21 

-83.  7 

5 

TOTAL 

1746618. 

inside 

site  ring  radars 

0  10000  1  10 

RANK 

NAME 

POPULATION 

CLOSEST  SITE 

DISTANCE 

AZIMUTH 

1 

NEW  YORK, NY 

10694632. 

RADAR2 

171. 79 

-157.  9 

2 

CHICAGO, 1 LL 

6220913. 

RADAR  5 

260. 72 

-75.  8 

3 

LOS  ANGELES, CAL 

6038771. 

RADAR3 

259. 56 

-84.  3 

4 

PHI LADELPHI A, PA 

4342897. 

RADAR8 

161. 32 

94.  7 

5 

DETROIT, MICH. 

3762360. 

RADAR5 

164. 86 

47.  1 

6 

SAN  FRANCISCO, CAL  2648762. 

RADAR9 

130.37 

-101. 9 

7 

BOSTON, MASS 

2595481. 

RADAR  2 

176. 36 

80.  8 

8 

PITTSBURGH, PA 

2405435. 

RADAR8 

259. 39 

-78.3 

9 

SAINT  LOUIS, MO 

2104669. 

RADAR  1 1 

227.94 

86.9 

10 

WASHINGTON, DC 

2076610. 

RADAR  8 

120. 50 

180.0 

10 

TOTAL 

42890464. 

cum  distribution  radars 

11  25 

RANK 

NAME 

POPULATION 

CLOSEST  SITE 

D 1  STANCE 

CUM 

NORMAL 

21 

SEATTLE, WASH 

1107213. 

RADAR1 

69.57 

1107213. 

0.0563 

12 

BALTIMORE, MD 

1803745. 

RADAR8 

84.02 

2910958. 

0.1481 

22 

KANSAS  CITY, MO 

1092545. 

RADAR11 

146. 97 

4003503. 

0. 2037 

20 

DALLAS , TEXAS 

1119410. 

RADAR4 

150.36 

5122913. 

0. 2606 

19 

PATTERSON, NJ 

1186873. 

RADAR2 

158. 36 

6309786. 

0.3210 

13 

NEWARK, NJ 

1689420. 

RADAR2 

172. 78 

7999206. 

0. 4069 

23 

SAN  DIEGO, CAL 

1033011. 

RADAR3 

201.69 

9032217. 

0.4595 

25 

INDIANAPOLIS, 1 MD 

944475. 

RADAR  5 

212. 43 

9976692. 

0.5075 

14 

MINNEAPOLIS, MINN 

1482030. 

RADAR10 

218.53 

11458722. 

0.5829 

18 

CINCINNATI, OH  10 

1268479. 

RADAR5 

249. 04 

12727201. 

0.6474 

11 

CLEVELAND, OHIO 

1909483. 

RADAR5 

252. 77 

14636684. 

0.7446 

17 

Ml LWAUKEF , Wl $ 

1278850. 

RADAR10 

258.94 

15915534. 

0. 8096 

15 

HOUSTON, TEXAS 

1418323. 

RADAR4 

265.37 

17333856. 

0.8818 

16 

BUFFALO, NY 

1306957. 

RADAR8 

356. 82 

18640800. 

0.9483 

24 

ATLANTA, GA 

1017188. 

RADAR6 

357. 79 

19657984. 

1.0000 

Fig.  VI-2.  DEANE  commands  for  Example  10. 
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VII.  CONVENIENCE  COMMANDS 


IN  PI  T  <p 

USE  -  This  command  allows  for  specifying  input  devices  other  than  the  default  device. 

It  is  useful  for  reading  large  or  frequently  used  files  of  DEANE  commands  from  stor¬ 
age.  Such  a  file  should  end  with  ’INPUT  5'  if  it  is  desired  to  return  control  to  the  de¬ 
fault  input  device. 

RESPONSE  —  (This  is  implementation  dependent  and  is  described  for  Lincoln  Labora¬ 
tory’s  implementation  of  IBM's  CP/CMS.  )  If  V  is  an  integer,  it  becomes  the  input 
device;  otherwise,  V  is  assumed  to  be  the  name  of  a  data  file  of  filetype  ’file',  and 
the  input  device  is  set  to  4  and  is  defined  to  correspond  to  file  V.  Reading  always 
starts  at  the  top  of  the  specified  input  file. 

ERROR  —  If  V  is  an  integer,  it  may  not  be  0,4,  6,  or  any  number  greater  than  12.  If 

V  is  a  filename  and  is  not  known  to  the  account,  an  error  message  is  printed  and  the 
command  is  ignored. 

OUTPUT  <p 

USE  —  This  provides  the  feature  for  writing  output  onto  a  data  file  so  that  the  file  off¬ 
line  can  then  be  printed. 

RESPONSE  —  (This  is  implementation  dependent  and  is  described  for  Lincoln  Labora¬ 
tory's  implementation  of  IBM’s  CP/CMS.  )  If  V  is  an  integer,  it  becomes  the  output 
device,  otherwise,  V’  is  assumed  to  be  the  name  of  a  data  file  of  filetype  ’file',  and 
the  output  device  is  set  to  33  and  is  defined  to  correspond  to  file  V.  Printing  always 
begins  at  the  end  of  the  existing  data  in  the  specified  file. 

ERROR  —  If  V  is  an  integer,  it  may  not  be  0,  4,  5,  or  any  number  greater  than  12.  If 

V  is  a  filename  and  is  not  known  to  the  account,  a  message  is  printed  to  this  effect, 
the  file  is  created,  and  execution  continues. 

LISTJSENSORS 

USE  —  This  command  affords  the  user  the  opportunity  to  obtain  a  listing  of  all  the  sen¬ 
sors  in  the  sensor  list,  together  with  their  positions,  boresights,  and  types. 

RESPONSE  —  The  listing  described  above  is  produced. 

ERROR  —  None. 

LISTJFARMS 

USE  —  This  command  allows  the  user  to  obtain  a  listing  of  all  the  farms  in  the  farm 
list,  together  with  their  positions,  populations,  and  types. 

RESPONSE  —  The  listing  described  above  is  produced. 

ERROR  —  None. 

DUMP  <p 

USE  —  The  two  commands  above  allow  the  user  to  view  the  contents  of  the  sensor  and 
farm  lists  during  execution.  This  command  permits  the  user  to  interrogate  the  contents 
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of  any  other  input  data  either  for  purposes  of  debugging,  review,  or  completeness  in 
the  output. 

RESPONSE  —  If  V  is  any  data  defining  command  other  than  'NEW_FARMS'  or  'NEW_ 
SENSORS’,  output  with  identification  will  be  produced  giving  the  values  associated  with 
that  command.  If  V  is  'ALL',  values  associated  with  every  data  defining  command 
will  be  printed  with  identification. 

ERROR  —  If  V  is  neither  ’ALL'  nor  a  DEANE  command,  a  message  will  be  printed  to 
that  effect. 

COMMENT  <p 

USE  —  This  command  allows  the  user  to  enter  a  line  in  the  output  for  purposes  of  iden¬ 
tification,  clarification,  or  documentation. 

RESPONSE  —  The  character  string  V  is  reproduced  on  the  output  device. 

ERROR  —  None. 

SPECIAL 

USE  —  This  command  provides  an  escape  into  the  code  for  the  DEANE  interpreter. 

RESPONSE  —  Control  passes  to  a  Fortran  subroutine  named  SPECIAL  which  the  user 
must  have  written,  compiled,  and  loaded  with  the  DEANE  interpreter.  As  the  inter¬ 
preter  is  written  in  Fortran,  all  the  data  and  subroutines  in  the  DEANE  interpreter 
are  accessible  via  this  command  if  the  user  understands  the  coding  of  the  interpreter. 

ERROR  —  None. 

STOP 

USE  —  This  command  provides  the  normal  exit  from  DEANE. 

RESPONSE  —  Control  returns  to  the  CMS  environment  from  DEANE. 

ERROR  -  None. 

Example  11 

All  the  data  which  can  be  specified  by  data  defining  commands  have  initial  values  at  the  be¬ 
ginning  of  a  session.  If  the  first  command  is  ’DUMP  ALL’,  as  seen  in  Fig.  VII-1,  these  initial 
values  are  printed.  In  printing  the  p  table,  *  ’s  result  because  the  initial  value  is  a  constant  p 

7 

of  10  which  is  too  large  to  fit  in  the  field  provided.  To  see  the  result  for  dumping  just  one 
command,  'DUMP  PARAMETER1  is  executed  as  seen  in  Fig.  VII-1. 
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Fig.  VII-1.  DEANE  commands  for  Example  11. 


VIII.  COMMAND  SUMMARY 


The  dependence  of  the  results  of  executing  a  command  on  the  results  of  executing  previous 
commands  is  seen  in  Table  VIII- 1.  The  vertical  column  contains  '  USEJSENSORS' ,  '  USE_FARMS' , 
’TRACKING  SENSORS’,  and  all  the  action  commands;  the  horizontal  entries  are  the  datadefining 
commands  and  'working  trajectory'  which  refers  to  the  data  defining  the  working  trajectory  as 
calculated  by  execution  of  either  ’KEPLER',  ’KEPLERP',  'EJECT',  or  'INTEGRATE'.  Reading 
horizontally  for  any  present  command,  an  'x'  indicates  that  the  present  command  does  depend 
on  the  data  defined  by  the  corresponding  horizontal  entry,  a  indicates  that  the  present  com¬ 
mand  may  or  may  not  depend  on  the  associated  horizontal  entry,  and  a  blank  indicates  that  there 
is  no  dependence. 

For  example,  executing  'INTERCEPT'  depends  on  'OUTPUT'  for  specification  of  where 
the  results  are  to  be  printed,  ' USE_JSENSORS'  for  specification  of  the  acquisition  sensors, 

'  FL YOUT_TABLE'  for  specification  of  interceptor  performance  characteristics,  'RESTRICT' 
for  system  parameters  limiting  intercepts,  'USE_F ARMS'  for  specification  of  an  interceptor 
deployment,  'TRACKING_SENSORS'  for  definition  of  the  tracking  radars,  and  the  'working 
trajectory'  for  defining  the  object  to  be  intercepted.  Data  from  '  RAI)AR_TYPE'  will  be  needed 
to  define  limiting  range,  elevation,  and  azimuth  for  acquisition  and  tracking  sensors  only  if 
some  radar  restriction  is  included,  i.  e.,  it  is  not  needed  if  there  are  no  sensors  in  use  and 
there  are  no  tracking  sensors  defined. 
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TABLE  VIII-1 

MATRIX  SHOWING  DEPENDENCE  OF  ONE  COMMAND 
ON  EXECUTION  OF  OTHER  COMMANDS 
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IX.  MATHEMATICS  AND  LOGIC 


A.  Basic  Coordinate  System 


It  proves  convenient  to  work  in  an  inertial  Cartesian  coordinate  system  whose  origin  is 
coincident  with  the  center  of  a  spherical  earth  of  radius  Re.  The  time  convention  adopted  is  that 
t  0  is  the  impact  time.  Time  is  referenced  as  time  before  impact,  and  it  has  a  negative  value. 
The  orientation  of  the  earth  relative  to  the  inertial  frame  is  defined  to  be  such  that  at  t  0  the 
x-axis  passes  through  0°  latitude  and  0°  longitude.  This  system  is  sketched  in  Fig.  IX- 1.  The 
position  of  a  point  at  latitude  <p ,  longitude  0,  and  height  h  above  the  earth  is  thus  given  by 


X 

cos  <p  cos  (0  +  cc  t) 
£ 

y 

=  (Re  +  h) 

cos  <z>  sin  (0  +  tc  t) 
g 

z 

sin  0 

where  a: 

g 


is  the  angular  rate  of  the  earth's  rotation  in  inertial  space. 


B.  KEPLER 

The  missile  is  assumed  to  move  in  an  inverse-square-law  central  gravitational  field  for 
its  entire  flight.^  That  is,  air  drag,  rocket  thrust,  lunar  and  solar  gravity,  asphericity  of  the 
earth's  gravitational  field,  and  the  earth's  orbital  motion  are  neglected.  Thus,  the  trajectory 
calculated  is  an  arc  of  an  ellipse  in  inertial  space. 

To  see  what  needs  to  be  known  to  specify  motion  on  a  trajectory,  refer  to  Fig.  IX-2.  The 
orientation  of  the  trajectory  plane  with  respect  to  the  inertial  system  is  defined  by  ft  (longitude 
of  the  ascending  node)  and  i  (inclination).  Motion  in  the  trajectory  plane  describes  an  ellipse 
and  the  spatial  relationship  is 

r  =  ~7 - - -  (1) 

m  1  —  e  cos  v 


t  The  analysis  in  this  subsection  was  performed  by  W.  C.  Kellogg. 
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|lH-t30»l 


Fig.  IX-2.  Keplerian  trajectory 
parameters. 


where  r  is  the  distance  from  the  center  of  the  earth,  p  is  the  semi-latus  rectum,  and  c  is 
m 

the  eccentricity.  The  angle  v  (true  anomaly)  is  measured  from  apogee  and  is  related  to  time 
via  Kepler's  equations.  The  angle  a;  (argument  of  apogee)  is  that  in  the  trajectory  plane  from 
the  ascending  node  to  apogee.  It  is  positive  if  apogee  is  in  the  northern  hemisphere,  and  neg¬ 
ative  if  apogee  is  in  the  southern  hemisphere.  With  these  definitions,  the  inertial  position  of  a 
projectile  is  given  by 


X 

m 

cos(o)  +  r)  cosfi  —  sin(o;  +  v)  sinfi  cosi 

y 

-  r 

cos  (to  +  v )  sin  £2  +  sin  (to  +  v)  cos  cosi 

J  m 

m 

sin  (to  +  v)  sini 

7  m 

L 

To  get  velocity  components  in  the  inertial  frame,  differentiation  of  Eq.  (1)  gives 


(2) 


ev  sin  v 


but  for  the  inverse-square-law  central  force  we  have 


r 


2 

m 


=  \pG 


where  G  is  the  universal  gravitational  constant  times  the  mass  of  the  earth.  Thus, 


m 

The  velocity  may  now  be  calculated  simply  as  the  time  derivative  of  Eq.  (2)  yielding 


r  = 

X 

m 

y 

r 

m 

X 

m 

y 

+  r  v 

-sin  (a;  +  v)  cosfi  -cos(u  +  y)  sin  £2  cosi 

-sin  (a;  +  v )  sin  ft  +  cos  (to  +  v)  cos  Q  cos  i 

m 

J  m 

z 

r 

m 

J  m 

z 

m 

cos  (to  +  v )  sini 

m 

m 

Thus,  to  define  a  Keplerian  trajectory,  we  need  to  find  values  for  the  constants  p, 
and  co  and  to  determine  the  functional  relationship  between  v  and  t. 


(3) 

(4) 


(5) 


c,  ft,  i, 
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1.  Values  of  Constants 


G  =  Universal  gravitational  constant  times  mass  of  earth 

=  0.3986094  x  106km3/sec2 
u;^  =  Angular  rate  of  earth  in  inertial  space 

=  0.72921  X  10  ^rad/sec 
He  =  Radius  of  spherical  earth 
-  6375.2km 


2.  The  Spherical  Triangle 

In  order  to  determine  the  plane  of  the  orbit  and  to  compute  the  distance  flown  by  the  missile, 
it  is  useful  to  consider  the  spherical  triangle  shown  in  Fig.  IX-3  with  vertices  at  the  launch 
point,  the  target,  and  the  North  Pole.t  The  sides  of  this  triangle  are  the  launch  point  and  target 
meridians  and  the  intersection  of  the  trajectory  plane  with  the  earth’s  surface. 


The  lengths  of  the  sides  meeting  at  the  North  Pole  are  the  complements  of  the  launch 
latitude  tp  j  and  the  target  latitude  and  are  known  from  input.  The  North  Pole  angle  a^p 
is  the  difference  between  the  target  east  longitude  0^  and  the  launch  point  east  longitude  0^  , 
plus  the  angle  through  which  the  earth  rotates  during  the  flight. 


a.  The  North  Pole  Angle  a^p 
The  angle  a^p  is  given  by: 


*NP 


^0T-0L 


cc  t,, 
g  F 


where  0T  and  0^  are  the  launch  and  target  longitudes,  a;  is  the  angular  rate  of  the  earth,  and 
Ij  I  g 

tp  is  the  time  of  flight.  The  difference  in  longitudes  is  used  as  an  initial  value,  and  the  cor¬ 
responding  time  of  flight  is  calculated.  This  time  of  flight  leads  to  a  new  value  of  a^p,  which 


t  The  analysis  here  relies  heavily  on  spherical  trigonometry.  The  same  results  can  be  de¬ 
rived  more  easily  using  vector  analysis  in  the  inertial  Cartesian  frame.  Since  the  subroutine 
which  performs  these  computations  is  based  on  the  trigonometric  analysis,  that  is  the  analysis 
discussed. 
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leads  to  a  new  time  of  flight,  etc.  Iteration  ceases  if  the  following  Cauchy-derived  criteria  are 
met:  let 


tj.  =  time  of  flight  calculated  from  this  iteration 
tj.  =  time  of  flight  used  in  this  iteration 
t^  =  time  of  flight  used  in  last  iteration 


(tp  tp  tj.’  are  three  consecutive  calculations  of  the  time  of  flight;  before  iteration  begins,  the 
time  of  flight  is  initialized  to  zero).  Then  we  have  converged  if 

(1)  the  last  two  times  are  nearly  equal:  1 1 ^  —  tj.  |  <  0.1  sec,  or 

(2)  the  last  two  times  are  somewhat  close,  and  solution  has  begun  to 
diverge:  1 1 ^  —  t |  <  0.5,  and  1 —  t |  >  1 1*.  -  tj.’  | . 


If  neither  (1)  nor  (2)  above  is  met  within  12  passes,  an  error  message  is  printed  and  the  cal¬ 
culation  is  aborted. 

In  the  fixed-earth  case  where  a;  =  0,  the  time  of  flight  has  no  influence  on  the  calculations. 
Thus,  only  one  pass  is  taken  and  iteration  ceases  after  the  first  set  of  calculations. 

The  angle  a^p  is  a  signed  quantity:  it  is  positive  if  the  target  at  impact  time  is  east  of  the 
position  of  the  launch  point  meridian  at  launch  time. 


b.  The  Flight  Path  Arc 
From  the  law  of  cosines  for  spherical  triangles, 

cos^p  =  sin<pj  sin  <p^  +  cos<pj  cos  (p^  cosa^p 


we  derive 


and 


(~f)  =  si"2  (“S— L)  +  COS<pL  COS<pT  sin2  ("IF) 

/"f\  .  2  (<p L  +  <PT\  ,  2  /aNP\ 

I  — I  =  sin  I - j - 1  +  cos<pj  cos  (prj,  cos  1-^ — I 


We  combine  these  two  equations  with 


,^F=  2  sin  (^)  cos(^) 


to  get 


This  is  superior  to  the  usual  method  of  calculating  sinr^, 

2  ,1/2 

sin  v  =  (1  -  cos  v ^  ) 
for  angles  of  close  to  90°. 


(6) 
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c.  Missile  Headings  Aj  at  Launch  and  Ap  at  Impact 

The  four-quadrant  angles  A^  and  Ap  specify  the  angles  from  north  to  the  projections  on 
the  earth's  surface  of  the  missile  velocity  vectors  at  launch  and  impact,  respectively.  If  the 
projection  points  east,  A  =  90°;  if  it  points  south,  A  =  180°;  if  it  points  west,  A  =  270°. 

From  the  law  of  sines,  we  have 


sin  Aj. 


NP 


sin 


cos  (pT 


Sm  aNP 

sinAT  =  IJmTTr  cosv 


I, 


and  from  the  law  of  cosines. 


cos  Aj.  = 


sin  ipr p  -  cos  sin  <p  ^ 

cos  (p  T  sin  v  _ 

L  F 


sin<pT  —  cosr  sin^rp 

a  I  j  r  1 

cos  A_  =  - . - 

T  cos^T 

Accordingly, 

sin  a^p  cos  <p  ^  cos  <pp 

Ap  tan  sin</?~  -  cosi^r  sin<pT 
1  r  L< 

sin  a  pcos^T  cos  o?T 

A  =  tan  1 _ — _ _  1 J 1 _ 

T  -sin^j  +  cos^p  sin  (p^ 


d.  The  Inclination  i 

The  ascending  node  is  the  node  at  which  the  trajectory  passes  from  the  southern  to  the 
northern  hemisphere.  The  direction  in  which  the  missile  moves  along  the  trajectory  is  specified 
by  the  inclination  1  from  the  equatorial  plane  to  the  trajectory  plane.  To  an  observer  in  space 
above  the  North  Pole,  0  ^  i  <  90°  for  a  counter-clockwise  trajectory,  and  90°  <  i 180°  for  a 
clockwise  trajectory. 

Since  i  is  an  angle  of  a  right  spherical  triangle, 

cosi  =  cos<pT  sin  A 
Lj  l 

and  substituting  for  sinA^  gives 


whence 


cos  i 


sina^pCos</?j  cos  <pp 
sinl,F 


i  =  tan 


_i  ,/sin2  -  (sin  a^.p  cos  (*>L  cos  (/>T)2 
sina^pCos<pj  cos  <pp 


since,  by  convention,  sin  i  ^  0. 
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e.  Longitude  of  the  Ascending  Node 


Since  no  sidereal  time  reference  is  given  in  this  program,  it  is  necessary  to  establish 
arbitrarily  that  the  longitude  of  the  node  will  be  given  as  earth  longitude  at  impact  time.  Thus, 
the  longitude  of  the  ascending  node  is  given  by 

0  =  ®T  “  aN  “  aNP 


where  0^  is  the  target  longitude,  a^p  is  the  converged  value  of  the  North  Pole  angle,  and  a^ 
(as  indicated  in  Fig.  IX- 3)  is  the  equatorial  arc  directed  from  the  ascending  node  to  the  foot  of 
the  launch  point  meridian. 

From  the  law  of  sines, 

sinAj  sin<pp 


and,  because  we  are  in  a  right  triangle. 


cos  a 


N 


cos  Aj 
sin  i 


so  that  substituting  for  tanAT  gives 

L 


N 


tan 


sin  a^p  sm<?L  cos^  cos  <pp 


sin  (pr 


cos  sin  (p j 


f.  Argument  of  Apogee  cu 

The  angle  lc  is  the  arc  along  the  trajectory  projection  from  the  ascending  node  to  apogee, 
positive  if  apogee  is  in  the  northern  hemisphere,  negative  if  it  is  in  the  southern  hemisphere. 

In  terms  of  quantities  in  Fig.  IX- 3,  it  is  given  by  co  -  -  v ^  because  if  apogee  occurs  between 

launch  and  impact,  v ^  <  0  by  convention.  The  angle  v  ^  comes  out  of  the  time-of- flight  iteration. 
The  angle  is  derived  from  standard  right-triangle  formulas 

sin  <p  j 


cosr^.  =  ctnAj  ctn  i 

Substituting  from  the  equations  for  cos  i  and  cosAT  gives 

1  j 

_a  sin  v F  sin  <p 

v  _  =  tan  — - - - 

N  sin  <p^  -  cos  cos  <p^ 


3.  Inplane  Trajectory  Parameters 

Figure  1X-4  illustrates  the  relevant  quantities  of  the  inplane  analysis,  from  which  we  are 
to  determine  the  parameters  of  the  ellipse,  position  of  apogee,  time  of  flight,  and  launch  and 
impact  velocity. 

Apogee  will  usually  occur  between  launch  and  impact;  but  it  need  not  if  launch  and  impact 
elevations  are  very  different,  as,  for  example,  if  the  missile  were  launched  against  a  surface 
target  from  a  satellite  in  orbit.  The  true  anomaly  r,  specifying  the  position  of  the  missile  in 
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its  trajectory,  is  measured  from  apogee  and  is  taken  positive  in  the  sense  in  which  the  missile 
moves.  The  direction  of  the  orbit  about  the  earth's  axis  is  specified  by  the  inclination  i. 

From  the  way  in  which  the  true  anomaly  v  is  defined,  we  have 


where  v  ^  is  the  true  anomaly  at  launch,  v  ^  is  the  true  anomaly  at  impact,  and  v ^  is  the  great 
circle  arc  traversed  by  the  missile  from  launch  to  impact. 

Inputs  to  the  inplane  analysis  are  r^  and  r^  (the  radii  at  launch  and  impact),  v p,  and  one 
other  datum.  From  these,  we  need  to  compute  the  semi-latus  rectum  p;  the  eccentricity  e; 
the  true  anomalies  v  and  v  at  launch  and  impact,  which  determine  the  position  of  apogee;  the 
time  of  flight  t^;  and  the  launch  and  impact  velocities  Vj  and  V^.  All  but  the  velocities  must 
be  determined  inside  the  iterative  loop  for  determining  tp  and  a^p. 

The  chordal  distance  c  from  the  launch  point  to  the  impact  point  is  an  important  parameter. 
The  basic  equation 


c 


2 


rL_2rLrT  COSl,F 


+  r 


2 

T 


may  be  altered  to  use  the  result  of  Eq.  (6)  for  added  accuracy  when  =  0 


<rT 


Trp) 


4rLrT 


sin 


Resides  the  quantities  v j  ,  and  v^t  the  fourth  datum  may  be  one  of  four  types: 

(a)  A  minimum  launch  velocity  (energy)  trajectory. 

(b)  0rp,  impact  velocity  angle,  measured  from  local  horizontal/^ 

(c)  |VT|,  magnitude  of  impact  velocity. 

(d)  h  ,  apogee  height. 


t  The  representation  of  the  angles  shown  in  Fig.  IX-4  has  them  measured  from  an  upward  di¬ 
rected  local  vertical,  and  this  convention  is  used  in  the  analysis  and  in  the  actual  computer  code. 
For  DEANE  input  and  output,  however,  the  angles  are  measured  from  the  local  horizontal.  Thus, 
Pj  (analysis)  =  90-0  (input  and  output)  and  0T  (analysis)  =  0T  (input  and  output)  +  90. 
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Before  developing  the  algorithms  for  these  four  types  of  solutions,  it  is  useful  to  have  an 
expression  for  the  velocity  in  the  trajectory  plane.  We  have 

V  =  rr  \  vv  v 


where  r  and  v  are  unit  vectors  in  the  radial  and  angular  directions,  respectively.  Substituting 
from  Eqs.  (1),  (3),  and  (4),  we  have 


(— e  sin  rr  +  (l-eeosr)  r] 


(8) 


a.  Minimum  Impact  Velocity  Solution 
If  we  take  the  magnitude  of  Eq.  (8),  substitute  v  =  v we  get 

V2  (1  -  2e  eosi'T  +  e2) 


and  minimize  it  (by  method  of  Lagrange  multipliers)  with  two  constraints  from  Eq.  (1)  at  launch 
and  impact,  respectively. 


p  -  rI  I1  -  e  eos(i'T  -  I'p)] 

(9) 

p  =  rT  (1  -  e  coskt] 

(10) 

We  arrive  at  the  basic  equation  for  the  minimum  energy  ease: 

2 

sin  r  (1  +  €  )  =  2e(sinr  -  sin  v  ) 

r  1  1  j 

(11) 

We  then  eliminate  p  from  Eqs.  (9)  and  (10)  to  get 

e(rT  cos  rT  -  r  cos  v^)  =  rT  -  r^ 

(12) 

Thus,  from  Eqs.  (11)  and  (12),  we  may  solve  for  r  in  terms  of  r  r^,  rT  ,  and  e: 

L  r  1  L 

r T  [(1  +  e2)  eosi'p-  2]  +  (1  -  e2) 

sin  y.  =  - =— - - - r-n - ; — ; -  sini' 

L  2e(rT  +  )  (1  -  cos  I'p)  F 

2  2 

rT  [(1  +  e  )  sin  k  -  2(1  -  cos  i>p)]  +  2rr  (1  -  cos  »>„) 

COS  v  T  - - — - - - - r-y-. - r - - - 

L  2€(rT  +  rL)  (1  -  cos^) 

2 

We  can  eliminate  from  the  two  equations  above,  leaving  a  quadratic  in  e  .  The  useful 
solution  of  the  quadratic  is 


3  —  cos  v  —  2(1  -  cos  v  ) 

£  r  r 


r  4-  r 

T  L 


1  +  COS  V 


If  =  180°,  the  original  quadratic  may  be  reduced  with  the  result 


r  —  r 
T  L 

r  4  r 
rT  L 
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If  c  <  0  (i^p  =  180°  and  >  r^),  then  we  set  -  0  and  apogee  is  placed  at  launch.  In  any  case, 
when  €  ~  0,  we  set  v  =  —v^/Z  and  the  trajectory  is  circular. 

Having  determined  e,  we  now  may  determine  v By  using  the  chordal  distance  c, 

-(c  +  r  cos  v  -  r  )  sin  v 

sin  v  =  - -n-T1 - 4 - L 

L  ec(l  +  cos  i>p) 


c  -  rT(l  -  cos 

COS  V  -  - 

L  ec 

From  Eq.  (8),  the  actual  minimum  velocities  are 

2  2G<C  “  rx  +  rL  COS  v 
1  T  "  rTrj  (1  +  cos  Vp) 

2  2G(C  “  rL  +  rT  cos^F) 

“  rxr  i  ^  +  cos  ^F) 

For  the  critical  cases,  180°  or  e  =  0,  the  velocities  are 


2GrT 


rT(rT 


+  rT 


2Grn 


rL(rT 


+  rT 


If  we  substitute  the  expression  for  cos^j  in  Eq.  (1),  we  get  an  expression  for  the  semi- 
latus  rectum 


rLrT 

p  =  — - —  (1  -  cos  v  F> 

The  angles  between  the  velocity  vector  and  the  upward  directed  verticals  at  launch  and 


impact  (/3j  and  /3rp)  are  determined  from  Eq.  (8).  For  any  angle  v , 

„  .  -1  e  cos  v  -  1 

8  =  tan  - : - 

M  e  sin  v 

At  launch,  we  substitute  for  cos^j.  and  sini^  in  Eq.  (13)  to  get: 

r  T  sin  v 

/3t  =  tan  — - - 

c  +  r cos^n  —  r_ 


(13) 


T  F  *  L 

Similarly,  for  impact,  we  use  Eq.  (11)  to  obtain  expressions  for  sin^^  and  cos^^,  and  get 

-1  rLsin,/F 

^T  —  ^  _  j,  CQSV  ^  p 

L  FT 
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b.  Flight  Path  Angle  at  Impact  Specified 
If  the  angle  is  specified,  we  can  obtain  from  Eq.  (8)  [see  Eq.  (13)] 


e  cos  v _  —  1 

tan/?-,  =  - : - 

T  €  sin  \> 


(14) 


which  can  be  solved  for  e  in  terms  of  /?T  and  the  unknown  angle  v  This  expression  for  e  may 
be  substituted  into  Eq.  (1)  to  yield 


p  =  r 


l(‘- 


cos  v 


L 


cos  vT  -  sin  tan/?T  ( 


and 


P 


COS  V  r 


cos  v T  -  sin  PT  tan/?T  , 


for  launch  and  impact,  respectively.  We  equate  these  two  expressions  for  p  and  use  Eq.  (7)  to 
get  an  equation  whose  only  unknown  is  v 

KrL  ~  rT^  tan^T  +  rL  sin  V sin  1  T  ”  r\Ji  “  COS  ^F*  COS  ^T  0 


=  tan 


-1 


rL(l  -  cos  ^F) 


(rL  -  rT)  tan/?T  +  r^  sin^ 


Equation  (14)  yields  e,  which  may  be  outside  the  interval  (0,  1)  if  an  improper  value  of  p  is 
input.  The  semi-latus  rectum  p  can  be  obtained  from  Eq.  (1).  p^  is  obtained  by  substituting 
v  j  in  Eq.  (13). 

c.  Impact  Velocity  Specified 

To  solve  the  case  in  which  |  |  is  specified,  we  first  substitute  from  Eq.  (1)  into  Eq.  (8) 

to  get 


V 


2  G( 6  -  2e  cos  +  1 ) 
T  ~  rT(l  -  €  cos  y  ) 


If  we  let  p  be  the  ratio  of  impact  velocity  to  escape  velocity. 


2  2 
V  r  V 
2  V  rp  A  X  T 


V  ‘ 


2G 


then,  solving  the  equation  for  for  e  cosi^  gives 


2  2 
€  +  (1  -  2  p  ) 

e  cos  v  =  - — - vf—1 

2(1  -P  ) 


(15) 
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Wc  then  substitute  v ^  v ^  in  Eq.  (12),  and  the  expression  for  e  cosr^  from  Eq.  (15)  in  that  . 

result,  which  yields 

[e2  +  (1  -  2pZ)]  (r„  -  rT  cosO  -  2(1  -  p2)  (r,_  -  r.  ) 

€  sin  =  - - - ± ^ - —  .  (If,) 

2(1  -  p  )  r i  sin  v 


We  can  now  combine  Eqs.  (15)  and  (16)  to  obtain  a  quadratic  in  e  of  the  form, 

2  4  2  2  2  1 

c  €  -  2Be  +  (D  +  B  )  — —  =  0 


which  has  solutions 


2  B  +  D 

e  =  - r - 


(17) 


and 


2  B  -  D 


(18) 


where  c  is  the  chordal  distance  from  launch  to  impact,  and 

B  =  2  [rL  sin  r^(l  -  p2)  ]2  -  2(rT  +  r ^  )  r^(l  -  cos  v^)  (1  -  p2)  +  c2 

2  2 
D  =  2(1  -  p  )rT(l  -  cos  v p,)  •  {rL(l  -  cos  [p  rL(l  +  cos  rp) 

+  2p(rT  -  rL  cos^p)  -  rL(l  -  cos^)]}  7 

The  angle  may  be  computed  from  Eqs.  (15)  and  (16),  and  the  semi-latus  rectum  p  from 

Eq.  (1).  The  velocity  at  launch,  VT  ,  is  determined  by  substituting  r  for  v  in  Eq.  (8).  Finally, 

L  L 

the  flight  path  angles  pT  and  (3 ^  may  be  determined  from  the  velocity  Eq.  (13). 

For  certain  values  of  Iv^1,  either  no  trajectory  exists,  or  more  than  one  exist.  The 
program  tests  for  all  possible  cases  as  follows. 

First,  a  minimum  energy  trajectory  is  solved,  and  the  minimum  velocity  V  is  com¬ 
puted.  If  |  |  <  V  no  trajectory  exists,  so  we  set  =  V^.^,  and  a  message  is  printed. 

Next,  p  is  computed,  and  tested.  If  p  ^  1,  the  given  velocity  is  greater  than  the  escape 
velocity  Vp.  The  trajectory  is  thus  hyperbolic  and  the  solution  is  not  within  the  scope  of  DEANE. 
So,  wc  set  =  0.99V£,  p  =  0.9801,  and  print  a  message  to  that  effect. 

For  values  of  between  and  V  ,  there  will  be  two  trajectories  [corresponding  to  the 

two  values  of  e2  given  in  Eqs.  (17)  and  (18)],  one  below  and  one  above  the  minimum  energy  tra¬ 
jectory.  The  one  below  is  chosen  [Eq.  (18)]  if  the  velocity  is  specified  with  a  minus  sign,  and 
the  one  above  [  Eq.  (17)]  is  chosen  if  the  velocity  is  specified  with  a  positive  sign. 


d.  Apogee  Height  ha  Specified 

The  radius  to  apogee  r.  is  given  by  r  .  =  h  +  Re.  At  apogee,  cos  v  -  1,  so  that 

ix  rx  a 


(19) 
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It  is  convenient  to  introduce  the  ratio 


Combining  Eq.  (19)  with  Eq.  (1)  at  v  -  v  yields 


€ 


1  ~  P 

1  -  P  COS  v I 


(20) 


When  substituted  in  Eq.  (1)  at  both  launch  and  impact,  this  value  of  e  yields  two  expressions 

for  p  in  terms  of  v  which  must  be  the  same  for  the  correct  value  of  v  ; 

Li  L 


P  = 


L 


(rA~rL)  cos^  1 

rA-rLCOS"L  J 


p  =  rT 


(rA~rL)  cosl,Tl 

rA  “  rL  COS"l  J 


Equating  these  and  substituting  v ^  ,  we  obtain 


rL  -  rT  =  (rL  -  prT  -  (1 -p)  rT  cosHpJ  cost^  +  -  p)  sini»F]  sini'k  .  (21) 

For  convenience,  the  quantities  in  square  brackets  in  Eq.  (21)  are  labeled  B  and  C,  respectively. 
Then,  we  define  angles  ip g  and  given  by 


sin  = 

sin  ip  ~ 

COS  Ip  = 

Then,  Eq.  (21)  becomes 
sirup  ^  = 

"l  = 

after  which  we  obtain  e 


AFT? 

B 

c 

J7  +  c2 

sin  (r  j  +  ip) 

tPs-<P 

from  Eq.  (20)  and  p  from  Eq.  (19). 


e.  Velocities  and  Time  of  Flight 

The  calculation  of  p,  c,  ,  and  is  carried  out  as  described  above  by  the  method  appro¬ 
priate  for  the  input  data.  Equation  (14)  and  an  analogous  one  for  give  the  impact  and  launch 
angles,  respectively.  An  equation  obtained  by  taking  the  magnitude  of  Eq.  (8)  provides  the  mag¬ 
nitudes  of  the  velocities: 

V  =  (l-2e  cos.-  +  e2)  . 
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To  calculate  the  time  of  flight,  it  is  convenient  to  introduce  a  new  angle  E  called  the 
eccentric  anomaly,  given  by 

E  / 1  +  e  ,  v 
tan  2  =  tan  2  • 


(22) 


The  eccentric  anomaly  is  zero  at  apogee  like  the  true  anomaly.  The  time  t  since  apogee  is 
given  by  Kepler's  equation  of  time: 

nt  =  E  +  e  sin  E  (2  3) 

where  n  is  the  mean  anomaly,  given  by 

n  = 


and  a  is  the  semi-major  axis: 


Equations  (22)  and  (2  3)  are  derived  in  standard  references  on  orbital  mechanics.  In  order  to 
obtain  the  total  time  of  flight,  we  evaluate  Eq.  (22)  for  v  -  v ^  and  v  -  in  order  to  get  values 
of  E  for  launch  and  impact;  compute  the  corresponding  times  t  and  t  from  Eq.  (2  3);  and  then 
subtract  ty  =  trp  -  tj  .  To  relate  the  true  anomaly  v  to  time  for  the  convention  that  t  =  0  at 
impact,  Eq.  (23)  is  modified  to  become 

n(t  +  t^)  =  E  +  e  sin  E  (24) 

from  which  E  may  be  calculated  by  a  Newton-Raphson  iteration  for  any  value  of  t.  The  true 
anomaly  then  follows  from  Eq.  (22)  as 

v  =  2  tan_1  (VrH tan  f )  •  (25) 


C.  EJECT 


The  problem  here  is  to  take  velocity  increments  in  the  coordinate  system  shown  in  Fig.  IX-5, 

add  them  to  the  velocity  vector  at  point  P  to  get  a  new  state  vector  in  the  inertial  coordinate 

system,  and  then  calculate  the  trajectory  parameters  from  this  new  state  vector. 

Let  (AV-,  AV-,  AV-)  be  the  input  velocity  increments  and  let  the  state  vector  at  P  (as 
x  y  z 

known  from  the  old  trajectory  via  Eqs.  (24),  (25),  (2),  and  (5)]  be  (x,  y,  z,  V  ,  V  V  ).  In  the 

x  y  z 


Fig.  IX-5.  Coordinate  system  for 
input  to  'EJECT'. 
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inertial  system,  the  vector  x  is  merely  the  velocity  vector  at  P  implying 

V_. 

1 

X  "  V 


where  V  =  /  V  ^ 
\f  x 

is  given  by 


2  Z  2  - 

+  Vy  +  Vz  .  The  vector  z  is  the  vector  normal  to  the  old  trajectory  plane  and 


-sini  sin  £2 
sini  cos  £2 
—cos  i 

Of  course,  y  =  z  x  x.  The  velocity  increment  in  the  inertial  system  is  then 


AV 

x 

AV 

_ 

y 

AV 

z 

’av_' 

x 

AV- 

y 

AV- 

z 

(26) 


giving  the  new  state  vector 

(x,  y,  z,  V  +  AV  ,  V  +  AV  ,  V  +  AV  )  e  (x,  y,  z,  x,  y,  z) 

x  x  y  y  ^ 

Straightforward  application  of  trajectory  relationships  now  gives 

.2 


-1 


G 


{l  -  jt  l(xy  -  yx)2  +  (xz  -  zx)2  +  (yz  -  zy)2]  a  1]1/2 


i  =  cos  1  {(xy  -  yx)  [(xy  —  yx)2  +  (xz  -  zx)2  +  (yz  —  zx)V  2) 


p  =  a(l  -  € 


(27) 

-1)1/2 

(28) 

•  i2il/2\ 

-  zx)  )  } 

(29) 

(30) 

where  r 


2,2,2 

+  y  +  z 


,  S  = 


.2  .2  *2 

x  +  y  +  z  ,  a  is  the  semi-major  axis,  €  is  the  eccentricity, 


i  is  the  inclination,  and  p  is  the  semi-latus  rectum. 

At  the  point  P,  we  can  calculate  v  from  Eq.  (1)  as 


,  /  r  -  p 
- 1  /  m 

r  e 
m 


Equations  (2)  and  (5)  give 

sin  (w  +  v)  =  z/sin  i 


cos  (u;  +  v  ) 


r  z  +  z 
m 


4 


e  sin  v 


\Gp 


sm  i 
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so  the  argument  of  apogee  a;  is  given  by 


-1  t  sin  (a;  +  v )  n 

cc  -  tan  [ - 7 - - — 7  3  —v 

cos (a;  +  v  )  1 

To  get  the  longitude  of  the  ascending  node  ft,  we  have  from  Eq.  (2) 


—  =  cos  (a,  +  v)  cos  ft  -  sin  (a;  +  v)  sin  ft  cos  i 


—  cos  (a;  +  v)  sin  ft  +  sin  (a;  +  v)  cos  ft  cos  i 
r 


which  is  two  equations  in  two  unknowns  —  cos  ft  and  sin  ft.  It  readily  follows  that  after  solving 
these  equations 


ft 


tan 


sin  ft  . 
cos  ft 


It  remains  to  fix  up  the  time  convention.  Using  Eq.  (1)  and  the  given  impact  altitude,  v ^ 

(the  true  anomaly  at  impact)  can  be  determined.  Similarly,  the  angle  v ^  at  ejection  can  be  cal¬ 
culated.  Equation  (22)  gives  the  associated  eccentric  anomalies  E^  and  E^.  It  follows  from 
Eq.  (24)  that 

E .  +  e  sin  E. 

A  1  1 


defines  the  time  from  apogee  to  impact  t^.  The  time  from  ejection  to  impact  is  thus 


E^  +  e  sin  E^ 


Let  tg  be  the  time  for  an  object  on  the  base  trajectory  to  travel  from  the  ejection  point  to 

impact.  The  difference  tp  —  tg  must  be  saved  and  taken  into  account  for  the  rotating  earth  case. 

Impact  is  still  defined  as  time  zero,  but  with  a  nonzero  time  difference  the  x-axis  no  longer 

points  through  0°  longitude  at  time  zero.  Thus,  if  impact  is  at  (Re  +  alt)  (u  ,  u  ,  u  ),  the  im- 

x  y  z 

pact  latitude  and  longitude  are  given  by 

-1 

impact  latitude  =  sin  (u  ) 


impact  longitude  =  tan 


Vf-V 


The  velocity  and  impact  angle  may  be  calculated,  knowing  the  trajectory  parameters,  from 
Eqs.  (8)  and  (13)  taken  at  impact. 


D.  R_I_AXIS 

This  problem  is  the  inverse  of  the  problem  for  EJECT.  Given  an  impact  latitude  (p ,  lon¬ 
gitude  O,  and  altitude  A,  and  given  an  ejection  time,  we  want  the  velocity  increments  which 
give  that  impact.  We  assume  the  earth  does  not  rotate  and  that  changes  in  the  trajectory  param¬ 
eters  are  small  so  that  linear  approximations  (via  Taylor  expansions)  are  applicable.  (As  things 
turn  out  with  the  linear  approximation,  there  are  a  set  of  velocity  increments  which  satisfy  the 
impact  conditions  and  the  answer  will  be  a  vector  instead  of  a  point.)  The  problem  is  solved  in 
the  inertial  coordinate  system  and  then  transformed  to  the  coordinate  system  of  Fig.  IX-5. 
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At  the  new  impact  point,  we  have 


V 

cos  cp  cos  0 

yl 

=  (Re  +  A) 

cos  cp  sinG 

_zl_ 

sin  cp 

For  the  base  trajectory  at  the  same  altitude,  we  can  get  from  Eq.  (1),  and  the  position  of 
that  point  in  the  inertial  frame  is 


^0 

y0 

zo 


=  (Re  +  A) 


cos  (cc  +  Vj^)  cos  SI  —  sin  (u  +  sin  S2  cos  i 


cos(w  +  sin  SI  +  sin(u  +  v  cos  SI  cos  i 


sin  (a;  +  v^)  sin  i 


(31) 


Defining  a  -  u:  +  r.  and  since  dr  =  0  at  the  impact  altitude,  differentiation  of  Eq.  (31)  gives 


0 

yo 


+  K 


da 

dS2 


(32) 


where 


zi 

zo 

■ 

di 

dx/da 

ax/0S2 

ax/ai 

B  = 

ay/a  a 

0y/8Sl 

ay/ai 

az/a  a 

0z/0S2 

8z/ai 

A 

What  we  need  now  is  the  matrix  which  expresses  the  relation  in  Eq.  (3  3)  below,  since  combining 
that  with  Eq.  (32)  will  specify  dV: 


(33) 


at  ejection 
v-v , 


da 

dV 

X 

dS2 

-  C 

dV 

y 

di 

dV 

at  impact 

z 

A 


0 


Toward  this  end,  consider  the  point  of  ejection  where  the  spatial  coordinates  are  identical 
for  the  base  and  new  trajectories.  Differentiation  of  Eq.  (29)  at  ejection  time  gives 


di  =  (8i/aV  di/dV  di/dV  ] 


dV 


dV 


dV 


(34) 
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Define  p  o;  +  rQ,  where  rQ  is  the  true  anomaly  of  the  base  trajectory  at  the  time  of  ejection. 
Differentiation  of  Eq.  (2)  gives 


dz  =  0 

dx  =  0 


d/3  =  —tan  p  ctn  i  di 


sin  p  cos  12  cos  i  +  sing  cos  p  sin 


cos  p  sin  sini  +  sinp  cos  p  cos  £1  sini  cos  i 


— — :) di 

■in  i  r>Ac  i  ' 


(35) 


The  relation  remaining  to  complete  Eq.  (32)  is  that  between  do  and  dp.  Differentiating  Eq.  (1) 

yields 

dp  +  r  cos  r  de 
er  sin  r 

and  differentiating  Eqs.  (27),  (28),  and  (30)  at  v  =  rQ  gives  dp  and  de  in  terms  of  dV  so 

*dV 


dr  =  [3r/3V  3r/3V  3r/3V  ] 

x  y  z 


dV 


dV 


Thus,  evaluation  of  dr  at  r  =  rQ  and  r  =  r^  allows 

dev  =  d(u;  +  r^)  =  do;  +  dr^  +  drQ  -  drQ  =  dp  +  dr^  -  drQ  .  (36) 

Equations  (34),  (35),  and  (36)  now  completely  specify  the  relation  in  Eq.  (3  3),  and  we  have 


(37) 


0 

Since  dfl  is  a  constant  times  di  (implying  the  matrix  BC  is  singular),  the  solution  to  Eq.  (37) 
for  dV  is  a  line  and  not  a  point.  We  seek  a  solution  to  Eq.  (37)  of  the  form 

dV  =  a  +  kb 

XX  X 

dV  =  a  +  kb 

y  y  y 


xl‘ 

xo 

"dV  ' 

X 

y* 

- 

^0 

=  BC 

dV 

y 

zi 

z0 

dV 

z 

dV  =  a  +  kb 
z  z  z 


(38) 


and  it  is  convenient  to  consider  Eq.  (39)  below  where  D  =  BC: 


d 

dV 

X 

X 

d 

=  D 

dV 

y 

y 

d 

dV 

z 

z 

(39) 
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By  setting  dV  =  0,  (d  ,  d  ,  d  )  equal  to  the  left-hand  side  of  Eq.  (37),  and  solving  Eq.  (39) 
z  x  y  z  _ _ 

yields  ax  and  a^  with  az  =  0  and  k  =  0  in  Eq.  (38).  Setting  (dx,  d^,  dz)  =  0,  dVz  =  1,  and  solving 

Eq.  (39)  gives  b  and  b  with  b  =  1  in  Eq.  (38)  completing  the  solution  for  the  vector  of  velocity 
x  y  z 

increments  in  the  inertial  frame.  Transforming  these  results  into  the  coordinate  system  of 
Fig.  IX-5  using  the  inverse  (transpose)  of  the  matrix  in  Eq.  (26)  gives 


AV— " 

X 

X 

X 

X 

y 

X 

z 

a  +  kb 

X  X 

a_ 

X 

lx 

1 _ ; _ 

AV  = 

AV- 

y 

= 

^x 

yz 

a  +  kb 

y  y 

= 

a_ 

y 

+  k 

b_ 

y 

AV- 

z 

z 

X 

z 

y 

z 

z 

k 

a_ 

z 

In 

_ i 

To  find  the  velocity  increment  with  minimum  magnitude 

AV2  =  (a-  +  kb-)2  +  (a-  +  kb-)2  +  (a-  +  kb-)2 

x  x  y  y  z  z 

setting  the  derivative  of  AV^  with  respect  to  k  equal  to  zero  gives 


k 


a_b_ 
x  x 


b2 

x 


+  a_b—  +  a_b. 

y> z  : 

+  b-  +  b 

y 


from  which  AV  may  be  evaluated. 


E.  INTEGRATE 

In  the  inertial  coordinate  system,  the  equations  of  motion  are 
x  =  -  G  -  j  ap  Vx 


y  = 


apVy 


z  =  — G  -  ~  apV z 
r 


(40) 


where  r  =  \/(x2  +  y2  +  z2),  V  =  */(?"  +  y^  +  'tF),  a  is  i/p.  p  is  the  atmospheric  density,  and  G 
is  as  before.  These  can  be  integrated  numerically  by  defining 


u  =  (x,  y,  z,  x,  y,  z) 


7  du  ,  ■  . 

f  =  dT  =  (x<  y-  z-  x'  y< 


so  that  the  differential  equation  (40)  becomes 
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f  = 


0 

0 

0 

_G_ 

3 

r 

0  - 


G_ 

3 


0 


0  -f  apV 


-  -  apV 


~Z  apV 


Using  the  predictor-corrector  scheme  of  Milne,  we  have  the  predictor 


u  .  +  \  At(8f  ,  -  4f  „  +  8f  .) 
p-4  3  p- 1  p-2  p-3 


and  the  corrector 


u  -  u  ,  +  4  At(  f  +  4f  .  +  f  J 
p  p- Z  3  p  p-1  p-2 


To  initialize  the  iteration  of  these  equations,  the  quantities  and  f  are  taken  at  the 
altitude  where  the  integration  begins  from  the  Keplerian  specification  of  the  working  trajectory. 
The  quantities  u_^,  f u_y  f  and  u_^  are  similarly  calculated  where  the  spatial  and  ve¬ 
locity  components  are  taken  from  Eqs.  (2)  and  (5),  and  the  accelerations  are  computed  from 
Eq.  (40)  with  a  -  0. 

The  results  for  the  state  vector  u  are  made  accessible  by  storing  them  in  a  table.  Ref¬ 
erences  to  the  integrated  part  of  the  working  trajectory  are  then  made  by  linear  interpolation  in 
this  table,  while  references  to  the  Keplerian  part  are  made  via  Eqs.  (2)  and  (5)  as  before. 

As  in  the  case  for  EJECT,  additional  consideration  must  be  given  to  the  time  convention. 
Consider  Fig.  IX-6  where  the  integration  begins  at  alt^  km  and  terminates  at  alt^  km.  The  time 
from  alt^  to  alt^  is  t^  sec  for  the  Keplerian  base  trajectory,  and  t.  sec  for  the  integrated  trajec¬ 
tory  with  |u  |  ^  |t^|.  With  a  rotating  earth,  this  time  delay  must  be  included  for  transforming 
longitude  into  the  Cartesian  frame  or  from  the  Cartesian  frame  into  longitude.  Also,  in  ref¬ 
erencing  the  Keplerian  part  of  the  trajectory  when  the  re-entry  portion  is  the  result  of  integra¬ 
tion,  |t.  -  t^  |  must  be  considered.  For  the  point  labeled  t^  in  Fig.  IX-6,  Eqs.  (2)  and  (5)  must 
be  evaluated  at  the  true  anomaly  v  corresponding  to  a  time  t  =  t^  +  |t.  ~  t^  |,  where  tQ  is  the 
time  before  impact  and  t^  <  0. 


Fig.  IX-6.  Trajectory  and  time 
change  from  'INTEGRATE'. 


F.  PRINT 

Given  a  sensor  at  latitude  <p ,  longitude  0,  and  height  h  above  the  earth,  and  given  a  time 
tQ,  the  problem  is  to  determine  the  quantities  listed  in  Table  IV-1. 
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1.  Time  Considerations 


Wc  have  to  be  concerned  about  time  in  two  systems  -  inertial  Cartesian  and  rotating  earth. 
Making  the  conditional  assignments  in  Eqs.  (41),  (42),  and  (43)  below,  we  have  the  position  of 
the  sensor  in  inertial  space  as  specified  in  Eq.  (44). 


Ati  = 


At2  = 


0 


ifto<0 


*0  ■ 


t^  —  time  of  flight  if  t^  >  0 

0  if  working  trajectory  not  from  'EJECT1 
flight  time  increase  due  to  'EJECT'  otherwise 
0  if  working  trajectory  not  from  'INTEGRATE' 
flight  time  increase  due  to  'INTEGRATE'  otherwise 

r  -  (Re  +  h)  u 


(41) 


(42) 


(43) 


where 


ug  -  (cos  <p  cos  (0  +  c^t^),  cos  ^  s*n  ^  +  Ct?  tjj),  sin  <p ] 


(44) 


t'  =  tn  +  At,  +  At0 
U  U  1  2 


Calculation  of  the  missile's  state  vector  is  similarly  conditionally  dependent  on  time  and 
the  nature  of  the  working  trajectory.  If  t^  >  — At^,  the  state  vector  is  found  by  linear  interpolation 
in  the  table  of  stored  integration  results.  Otherwise,  it  is  taken  from  Eqs.  (2)  and  (5).  (As  noted 
at  the  end  of  Sec.  IX-E,  if  the  re-entry  portion  is  obtained  from  'INTEGRATE',  Eqs.  (2)  and  (5) 
are  evaluated  at  t^  +  At^  not  at  t^  for  t^  <  — At^.)  At  any  rate,  we  get  a  state  vector  Eq.  (45)  in 
the  inertial  system,  and  the  problem  is  now  one  of  combining  Eqs.  (44)  and  (45): 


(x,  y,  z) 


r  =  (x,  y,  z ) 
i  m  *  ' 


(45) 


2.  Range,  Azimuth,  and  Elevation 
The  radar  range  is  simply  given  by 


It  is  useful  to  define  a  unit  vector  u^  along  the  line  of  sight,  and  also  two  other  unit  vectors  u^ 
and  u^  which  point  east  and  north  from  the  sensor  and  are  in  the  plane  tangent  to  a  spherical 
earth  (of  radius  Re  +  h)  at  the  sensor  site.  The  vectors  ug,  u^,,  and  u^  then  define  a  right- 
handed  coordinate  system.  These  definitions  are: 
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(46) 


u  = 
r 


r  -  r 
m  s 


u E  [—  sin  (0  +  cos  (0  4  u gt q ) ,  0] 

=  |-cos  (0  4  o^t^)  s*n  <f>f  —3111(0  4  sincp,  cos<p] 

The  azimuth  a  from  north  and  the  elevation  e  are  now  given  by 


.  . / u  • 

,  - 1  .  sin  a  .  .  -1 /  r  E 

a  =  tan  (  )  =  tan  1 

u  •  u,T 
r  Pm 


e  =  sin  1  ( ur  •  ug) 


3.  Time  Derivatives  of  Range,  Azimuth,  and  Elevation 


It  must  be  remembered  that  ug,  u^  ,  and  u^T  rotate  with  the  earth.  To  obtain  the  velocity 
seen  by  an  observer  in  the  rotating  coordinate  system,  we  use  the  operator  relation 


_d_ 

dt 


rotating 


d_ 

dt 


fixed 


(47) 


which  is  derived  in  texts  of  rigid-body  mechanics.  In  Eq.  (47),  u  is  (0,  0,  a^).  The  missile 
velocity  v  as  observed  by  the  sensor  is  then  given  by 


v=r  —  a;  X  r  =(x  4  a;  y  ,  y  ~  a;  x  ,  z  ) 

o  o  m  g J  m  J  m  g  m  m 


(48) 


In  order  to  obtain  time  derivatives  of  scalar  range,  azimuth,  and  elevation,  two  more  unit 
vectors  are  introduced,  u^  is  a  unit  vector  in  the  horizontal  plane  pointing  in  the  direction  of 
increasing  a,  and  ug  points  in  the  direction  of  increasing  e  and  is  normal  to  the  unit  vector  u^ 
of  Eq.  (46)  in  the  plane  containing  u^  and  the  unit  vector  ug: 

Ua  =  UF  COS  a  ~  UN  s*n  a 

u  -  —  u^  sin  a  sin  e  —  uAT  cos  a  sin  e  4-  u  cos  e 
e  E  N  s 

where  Up,  u^,  and  u^  are  the  unit  vectors  at  the  sensor  site  previously  defined. 

In  terms  of  this  set,  we  have 


v  ~  r  u  4  r  a  cos  e  u  4  r  e  u 
o  r  o  a  o  e 


(49) 


where  rQ  is  the  Doppler  velocity  observed  by  the  sensor.  Thus,  the  equations  needed  to  convert 
v  of  Eq.  (49)  into  time  derivatives  of  radar  observables  are 

r  =  v  •  u 
o  r 


r  cos  e 
o 


(50) 
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4.  Range  Acceleration 


The  easiest  way  to  get  the  range  acceleration  r*o  is  to  differentiate  Eq.  (50),  giving 


r  =  v  •  u  +  v  •  u 
o  r  r 

Differentiation  of  Eq.  (48)  by  application  of  Eq.  (47)  gives 

v  =  [x  +  cu  (2  y  -  u ;  x  ),  y  -  cc  (2x  +  cc  y  ),  z  ] 

m  g  J  m  g  m  J  m  g  m  gJm  '  m 

where  (x^,  y  ,  zm)  is  defined  in  Eq.  (40),  and  differentiation  of  Eq.  (46)  gives 


u  = 
r 


(x  +  o;y,y-a;x,z)  u  ,  .  , 

—  g*  s'  -  g  s'  r  I—  —  ,  1 

— - — -  —  — ; — t  I  r  •  r  +  a;  (x  y  — y  x) 

I^T  |  |7  l  °  m  g  mJ  s  ■  m  s  J 


5.  Altitude  and  Ground  Range 

The  altitude  above  the  earth  is  simply 

alt  =  I  r  I  —  Re 
1  m 1 

The  ground  range  at  time  t  is  defined  as  the  distance  measured  along  the  earth's  surface 
from  the  sensor  to  a  point  on  the  earth  directly  beneath  the  missile,  so 

ground  range  =  Re  cos 

6.  Velocity  and  Body  Aspect  Angles 

The  angle  between  the  line  of  sight  and  the  missile's  velocity  vector  is  simply 

velocity  aspect  angle  =  cos 

The  body  aspect  angle,  however,  depends  on  whether  orientation  is  specified  via  'POINTING' 


or  'TUMBLE'. 


a.  Orientation  Specified  via  'TUMBLE' 

If  the  tumbling  rate  is  and  we  assume  random  initiation  of  the  tumbling,  all  we  need 
is  a  vector  which  traces  out  a  circle  in  the  trajectory  plane.  This  is  easily  defined  in  the  in¬ 
ertial  coordinate  system  by  using  Eq.  (2)  with  cc  +  v  replaced  by  cu^t.  Thus, 

p^  =  (coso),pt  cos  ft  -  sinw^t  sin  ft  cosi, 

cos  o;Tt  sin  ft  +  sina;Tt  cos  ft  cosi,  sino^t  sini) 


and 


-1  — 

body  aspect  angle  =  cos  ( p 


u  ) 
r 
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b.  Orientation  Specified  via  'POINTING 


The  altitude,  latitude,  and  longitude  specified  in  'POINTING'  determine  a  vector  a  in  the 
inertial  coordinate  system.  This  vector  is  simply  found  by  turning  the  altitude  into  time,  find¬ 
ing  the  coordinates  of  the  missile  rm  and  the  coordinates  of  the  latitude  and  longitude  r^  at  this 
time  so  that 


a 


(<*,  0,  y) 


Referring  to  Fig.  IX-7  we  need  to  get  the  vector  b  which  is  due  to  spinning  at  an  angle  ip  from 
a  and  with  a  rate  te  .  This  determines  c  -  a  +  b  so  that  the  aspect  angle  is 


body  aspect  angle  -  cos 


Fig.  IX-7.  Coordinate  systems  used  for 
finding  'PRINT'  mnemonic  'POINTING'. 


To  express  b,  define  the  right-handed  orthogonal  coordinate  system  (x',  y',  z')  such  that 


z '  =  (a,  0,  y )  e  a 

y1  =  (-/?,  Of,  oWflf2  +  02 

x'  =  y'  x  z'  -  [  ory,  py,  —  (a2  +  P^)]/J  a2  +  /J2 
The  length  of  b  is  tan 4,  and  it  rotates  in  the  (x',y')  frame  such  that 


(51) 


b  ,  ,  ,  =  (tani/j  cos  or  t,  tanip  since  t,  0) 

x  3  J  >  z  s  s 


Transforming  b  into  inertial  frame  coordinates  using  Eq.  (51),  we  have 


x,  y,  z 


7  Q2  +  /32 


- 1 

p 

1 

tan^  coscegt 

0y  a 

tan^  since  t 

,  Z  2, 

s 

-(a  +  0  )  0 

60 


7.  Cross  Section  and  S/N 

The  cross  section  is  found  by  interpolation  in  the  cross-section  table  specified  by 


1  RV_TYPE'  using  the  angle  determined  as  the  body  aspect  angle. 
With  this  cross-section  a,  the  S/N  is  evaluated  as 


where  E,  T,  O,  f  and  L  are  defined  via  * RADAR_TYPE',  and  R  is  the  range  in  kilometers. 


8.  Latitude  and  Longitude 


To  find  the  latitude  and  longitude  directly  beneath  the  missile,  consideration  must  be  given 
to  time  as  discussed  in  Sec.  IX-F-1.  Taking  the  inverse  of  the  relation  between  a  point  on  the 
rotating  earth  and  the  inertial  frame  which  is  expressed  in  Eq.  (41),  we  have 


latitude  =  sin 


where  t^  is  defined  in  Eq.  (44). 

G.  VISIBLE 

We  are  given  a  working  trajectory  and  a  sensor  at  latitude  <p ,  longitude  0,  height  h  with 
an  azimuth  borcsight  6.  It  is  not  easy  to  find  the  times  at  which  a  projectile  on  the  trajectory 
is  first  visible  and  is  last  visible  as  limited  by  range  R,  elevation  E,  and  azimuth  scan  ±  A  in 
an  analytic  fashion,  but  it  can  be  done.  A  search  is  certainly  less  aesthetic  than  application  of 
three  numerical  nonlinear  equation  solvers,  but  it  proves  to  be  speedier  in  terms  of  execution 
time.  Accordingly,  a  search  is  used  to  perform  the  calculation,  with  some  analysis  applied  to 
determine  starting  points. 

We  wish  to  find  two  times  at  which  to  start  searches  for  visibility  —  one  time  for  the  early 
visibility,  and  one  for  the  late  visibility.  To  do  this,  we  assume  that  the  whole  trajectory  is 
Keplerian  and  work  with  the  true  anomaly  v  as  an  equivalent  to  time.  Initially,  the  two  points 
are  v ^  and  corresponding  to  the  launch  and  impact.  Assuming  that  E  =  0,  we  can  find  the 
two  points  where  the  plane  specified  by  the  sensor  position  and  E  =  0  intersects  the  trajectory, 
implying  associated  values  of  v  ^  and  v ^  with  v ^  A  similar  computation  on  the  range  re¬ 

striction  leads  to  the  problem  of  finding  (or  approximating)  the  intersections  of  an  ellipse  and 
a  sphere.  This  implies  values  of  and  v ^  with  <  v ^  To  start  the  search  for  first  visibility, 

the  time  corresponding  to  v  =  max(^  ,  Vy  v  is  chosen.  Similarly,  to  start  the  search  for 
last  visibility,  the  time  corresponding  to  v  min(^fp,  v ^)  is  chosen. 


tin  degenerate  cases,  i.e.,  one  or  no  intersections,  the  two  anomalies  are  set  to  v  and 
respectively. 
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1.  Analysis  on  Elevation 

The  plane  of  0°  elevation  is  given  by 

x  y  z 

o  ,  ^o  ,  o 

- 2  x  +  - Z  y  +  - 7 

(Re  +  h)6  (Re  +  h)  (Re  +  h) 


z  =  1  =  nx  +  /3y  +  yz 


(52) 


where 


(x  y  f  zq)  =  (Re  +  h)  (cos  <p  cos  0,  cos  <p  sinG,  sin  <p ) 


Introducing  a  Cartesian  coordinate  system  in  the  trajectory  plane  such  that 


4  =  r  cos  v 
m 


ri  =  r  sin  v 
m 


we  have  the  relationship 


X 

V 

y 

=  T 

V 

z 

inertial 

where 


T  = 


cos  ce  cos  S3  -  sin  co  sin  S3  cos  i  -sin  co  cos  S3  -  cos  co  sin  S3  cos  i 
sin  S3  cos  co  +  since  cos  S3  cos  i  —since  sin  S3  +  cos  ce  cos  S3  cosi 


since  sini 

Substituting  in  Eq.  (52)  gives 

(a,  P,  7) 


X 

4 

y 

=  (a,  p,  y)  T 

V 

z 

cos  ce  sin  l 


=  1 


where 


M  + 


A  oTjj  +  0T21  +  yT31 


*=  aT2 1  +^T22  +  ?T32 


^  p  cos  r 


€  COS  y 


=  P  sin  k 
1  —  £  COS 


(53) 
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To  solve  Eq.  (5  3)  for  v,  we  substitute  sin 
tion  is 


in  v  -  \l l  - 2 


cos  v  and  get  a  quadratic  whose  solu- 


„„„  ,,  _  e  +  Ap  ±  <pp  n/e  2  +  A2p2  +  2 Ape  +  i^2p  -  1 

- 2 - 2  2  - 2~ 2 - 

e6  +  A  p  +  2 Ape  +  ip  p 


(54) 


From  Eq,  (53), 


sin  y  =  +_APL?.°i^ 

sip 


so  that 


r  =  tan 


sin  v 
cos  v 


Since  Eq.  (54)  gives  two  values  of  cos^,  two  values  for  r  are  obtained. 


2.  Analysis  on  Range 

The  intersection  of  a  sphere  of  radius  R  centered  at  the  sensor  and  the  trajectory  plane 
is  a  circle.  To  describe  this  circle,  the  equation  of  the  trajectory  plane  is 

Ax  +  By+Cz  1 
O  or  O 


where 


(A  ,  B  ,  C  )  =  (— sin  ft  sini,  cos  ft  sin  i,  — cosi) 
o  o  o 


With  (xo,  y  ,  zq)  given  by  Eq.  (52),  the  distance  from  the  sensor  to  the  trajectory  plane  is 

d  =  |A  x  +  B  y  +Cz  I 
1  o  o  oJ  o  o  o' 


so  the  radius  of  the  circle  is 

*  g 


=  Vr2  -  d2 


The  center  of  the  circle  is  at  (£  q,  tjq)  as  expressed  in  Eq.  (55)  below 

k  =  Ax  -By  -  C  z 
o  o  o^o  o  o 

x.  kA  +  x 
1  o  o 


y  .  =  kB  +  y 
J 1  oJo 


z  .  =  kC  +  z 
1  o  o 


V 

cos  u;  cos  ft  cos  a;  sin  ft 

—sin  ic  cos  ft  —sin  a;  sin  ft 


sin  a; 
sin  i 

cos  a; 


sin  i 


(55) 
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V 


Fig.  IX -8.  Geometry  for  bounding  time 
from  range  limits  in  'VISIBLE'. 


f 


In  the  trajectory  plane,  we  have  the  situation  shown  in  Fig.  IX-8.  Being  satisfied  to  only 
bound  the  angles,  we  have  Eq.  (56)  which  gives  values  for  and  v  ; 


f  = 


/7T“  2 

N  O  O 


"o  =  tan’1(7,o/f-  ^o/f) 


v~  -  v  -  sin  1  (-£•) 
3  o  f 


-  \  a 

vA  =  v  +  sin  (■#-) 
4  o  f 


(56) 


3.  The  Search 

It  is  assumed  that  the  period  of  visibility  is  continuous  in  time.  The  idea  of  the  search  is 
then  to: 

(a)  Using  t  =  function  (^  )  as  a  starting  point,  find  a  time  t  when  the 
missileeis  visible  by  using  a  gross  search. 

(b)  Push  t  as  far  backward  in  time  as  possible  with  a  fine  search  to 
get  thee earliest  visibility. 

(c)  Using  t^  =  function  (^)  as  a  starting  point,  find  a  time  t^  when  the 
missile  is  visible  by  using  a  gross  search. 

(d)  Push  t^  as  far  forward  in  time  as  possible  with  a  fine  search  to  get 
the  latest  visibility. 

The  logic  for  steps  (a)  and  (b)  above  is  shown  in  Fig.  IX-9,  and  the  logic  for  steps  (c)  and  (d) 
is  similar. 

H.  INTERCEPT 

Here  the  basic  problem  is  one  of  defining  a  fast,  complete,  logically  complex  algorithm. 

The  structure  of  the  algorithm  is  shown  in  Fig.  IX- 10.  The  actual  computer  code  is  slightly  more 
subtle  because  of  logic  to  speed  execution,  but  it  is  basically  that  seen  in  Fig.  IX- 10. 

I.  CONTOUR 

As  for  'INTERCEPT',  the  problem  here  is  to  define  a  fast,  complete,  and  logically  complex 
algorithm.  The  basic  structure  of  the  algorithm  is  shown  in  Fig.  IX-11,  although  the  actual 
code  is  much  more  complicated  in  an  attempt  to  speed  execution. 


64 


Similar  search  for  t| 
using  initial  t. 

from  ly  |l8-3  -I307s| 


Fig.  IX- 9.  Search  logic  for  ’VISIBLE’. 
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Fig.  DC-10. 


Basic  logic  for  ’INTERCEPT'. 


66 


Fig.  IX-11.  Basic  logic  for  'CONTOUR'. 
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Note  in  Fig.  IX-11  that  many  trajectories  are  calculated.  If  the  contour  is  for  vacuum 
trajectories,  each  trajectory  is  calculated  using  the  mathematics  for  'KEPLER’.  However,  if 
the  contour  is  to  be  generated  including  atmospheric  slowdown,  an  integrated  re-entry  portion 
is  calculated  once,  transformed,  and  then  stored  in  a  table.  These  data  are  then  transformed 
to  correspond  to  different  impact  points  (for  different  trajectories),  and  a  Keplerian  exoatmos- 
pheric  portion  is  calculated  and  tagged  on  to  the  re-entry  portion.  To  see  the  transformations, 
assume  the  earth  is  nonrotating  and  let 


launch  point  =  PQ  =  (xq,  yQ,  zq) 


actual  impact  point  =  p^  =  (x^,y^,z^) 


for  the  integrated  trajectory  calculated  using  the  mathematics  for  'INTEGRATE'.  To  store 
these  data,  define 


P. 


x 


P  =  P  x  P 
*y  z  x 


The  trajectory  plane  is  the  p  ,  p  plane,  with  p  pointing  through  the  impact  point. 

x  y  x 

Performing  the  transformation  (57)  below  on  (x,  y,  z)  and  (V  ,  V  ,  V  )  gives  a 

x  y  z 

P  P  P 
xx  rxy  xz 


a 


Pyx  Pyy  Pyz  ainertial 


(57) 


P 


four-element  (p  ,  p  ,  V  ,  V  )  table  to  store.  To  make  these  data  correspond  to  a  different 
_  Ky*  px*  py 

impact  point  p^  =  define 


q. 


Lz 


poXp 
^  l  *o 

X  p 


q  =  q  X  q 
My  Mz  Mx 
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For  each  point  stored,  we  perform  Eq.  (58)  below  to  get  (x,  y,  z)  from  (p^,  p^,  0)  and  (V^,  V  ,  Vz> 
from  (V  V  '  0). 

x  y 


qxx 

qyx 

^ZX 

°  inertial 

q 

^xy 

qyy 

qzy 

O' 

p 

(58) 

qxz 

qyZ 

qzz 

In  Fig.  IX-11,  there  are  several  references  to  getting  flyout  times.  There  is  some  logic 
coded  here,  not  shown  in  the  figure,  which  terminates  searches  for  intercepts.  If,  in  two  suc¬ 
cessive  calls  for  flyout  times,  say  t^  and  t^  (with  t^  <  t^),  neither  intercept  is  acceptable  (be¬ 
cause  of  reach  or  flyout  time)  and  the  range  to  the  trajectory  from  each  farm  in  use  is  increasing 
from  time  t^  to  time  t^,  then  no  acceptable  intercepts  will  be  found  by  searching  further  back 
in  time  and  the  search  is  terminated. 
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